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The absorption spectra for twenty-five compounds, of 
which seventeen are typical of the important components 
of gasoline and which include both chain and ring mole- 
cules, all in the vapor state, have been examined through- 
out the spectral range 1y to 30u, by using prisms of NaCl, 
KCl and KBr. The percentage absorption is shown 
graphically in the charts which follow as function of the 
wave-length for specified amounts of material. The curves 
indicate certain distinctive features characteristic of each 
different molecule, as well as bands common to several or 


all the compounds. The resolution was not sufficient to 
separate rotational lines in any of the bands. In addition to 
the hydrocarbons, eight other substances, tetraethyl lead 
(the only material studied in the liquid state), tetraethyl 
germanium, tetramethyl tin, and tetramethyl lead and also 
aniline, dicyclopentadiene, methyl nitrate and diethyl 
peroxide were investigated. They were included in the list 
on account of their influence upon the detonation of fuel in 
the gasoline engine. 


INTRODUCTION 


HE following paper presents the results of 

a series of absorption measurements which 
were carried out under the auspices of the 
Department of Engineering Research at the 
University of Michigan for the Research Labora- 
tory of the General Motors Corporation. The 
materials examined include the types of hydro- 
carbons which make up the principal components 
of gasoline, and certain other materials the 
interest of which lies in their effect on the 
detonation or “knock”’ of fuels.in the gasoline 
engine. In Table I the numbers refer to the 


TABLE I. 
Chart Chart 
Material number Material number 
Cyclohexane 1 Aniline 13 
Diiso butylene 2 Tetraethyl germanium 14 
Normal hexane 3 Tetramethy! tin 15 
2, 3 Dimethyl! butane 4 Tetramethy! lead 16 
Cyclohexene 5 Benzene 17 
1 Hexene 6 Toluene 18 
1 Octene 7 3 Heptene 19 
Methyl cyclohexane 8& 2 Methyl pentane 20 
Normal pentane 9 2, 2, 4-Trimethy! pentane 21 
Dicyclopentadiene 10 3 Methyl pentane 22 
Gamma heptine 11 2, 2 Dimethylbutane 23 
Diethyl peroxide 12 Tetraethy! lead 24 
Methyl nitrate 25 


charts giving the absorption data in graphical 
form. 


APPARATUS 


Two different spectrometers were used in 
mapping the following spectra. The first was of 
the Wadsworth type, equipped with a 60° NaCl 
prism, mirrors 12 cm in diameter and of 50 cm 
focal length, and a thermopile constructed by 
Firestone.' The small thermopile case was evacu- 
ated by a mechanical pump and activated 
charcoal, and was joined to a Moll galvanometer 
of d’Arsonval type. The combination was sensi- 
tive enough so that, with slits as small as 0.1 mm, 
it was possible to get good deflections out to 
4u, and with slits at 0.5 mm out to 14y, 
without a relay. At 6y the slits were usually 
0.2 mm wide, and this corresponds to a spectral 
range of 0.074. Deflections were observed on a 
translucent scale. After the angle of the prism 
had been measured, and the position was deter- 
mined in which light, reflected from the first 


1F. A. Firestone, Radiation Thermopile Design, Rev. Sci. 
Inst. 1, 630 (1930). 
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face of the prism, returned to the entering slit, 
it was possible to compute the figures for a 
calibration curve from the indices of refraction 
for NaCl? 

For part of the work the recording spectro- 
graph described by Randall and Strong* was 
available. This instrument was equipped first 
with a 30° prism of potassium chloride with a 
surface silvered by evaporation. With this 
arrangement it was possible to follow the spectra 
to 20u. It was used instead of the Wadsworth 
spectrometer for some of the materials (as 
indicated on the charts), and nearly all the 
spectra were extended by it beyond the range of 
the rocksalt prism. The measurements were 
carried still further by the use of a potassium 
bromide prism of 60° in the same apparatus, 
then equipped with a stationary and a moving 
mirror. The indices of refraction of KBr are 
now available,* but the recording spectrometer 
was calibrated empirically by means of absorp- 
tion bands of known wave-lengths, principally 
those of water vapor, carbon dioxide, ammonia, 
acetylene and the methyl halides. These wave- 
lengths have been determined in this laboratory 
by the use of gratings.® The absorption band of 
acetylene of 13.54 was photographed on nearly 
every record in order to have a check on the 
calibration. 

The absorption cell used with the Wadsworth 
spectrometer was provided with windows of 
NaCl for most of the work, but mica was used 


? These numbers, as determined by Paschen and others, 
are now accessible, with references to original sources, in 
Das Ultrarote Spektrum, Cl. Schaefer and F. Matossi, 
Springer, Berlin, pp. 45-47 (1930). 

7H. M. Randall and John Strong, A Self-Recording 
Spectrometer, Rev. Sci. Inst. 2, 585 (1931). 

*E. Gundelach, Die Dispersion von KBr Kristallen im 
Ultraroten, Zeits. f. Physik 66, 775 (1930). 

5E. K. Plyler and W. W. Sleator, Further Study of the 
Absorption of Infrared Radiation by Water Vapor, Phys. 
Rev. 37, 1493 (1931). 

P. E. Martin and E. F. Barker, The Infrared Absorption 
Spectrum of Carbon Dioxide. Phys. Rev. 41, 291 (1932). 

E. F. Barker, Molecular Spectrum of Ammonia; II, 
Double Band at 10u. Phys. Rev. 33, 684 (1929). 

Charles F. Meyer and Aaron A. Levin, The Infrared 
Absorption Spectra of Acetylene, Ethylene and Ethane. Proc. 
Nat. Acad. Sci. 13, 298 (1927). 

Willard H. Bennett and Charles F. Meyer, Infrared 
Absorption Spectra of the Methyl Halides. Phys. Rev. 32, 
888 (1928). 


where suitable, and some spectra were taken 
with windows of cellulose, made by dropping a 
little lacquer on water and lifting off the resulting 
film. These windows, however, were very fragile 
and proved pervious to the lighter hydrocarbons. 

With the self-recording spectrometer the ab- 
sorption cell, which had bromide windows 
capable of withstanding pressure, was inside the 
exhausted case and was pumped out separately 
through a small metal tube. It was connected 
to a mercury manometer so that the pressure 
of the inclosed gas could be read directly. In 
most cases records were taken with several 
different pressures, as best brought out the 
details of the absorption. A flask containing the 
liquid hydrocarbon was sealed to a second metal 
tube leading to the absorption cell; and opening 
a stopcock permitted more or less evaporation 
into the evacuated cell. 

With the Wadsworth spectrometer the ab- 
sorption cell was not pumped out, but some of 
the liquid to be studied was exposed inside the 
cell so as nearly to saturate the space within. 
The absorbing effect of constituents of the air 
was eliminated by moving the cell into and out 
of the beam of radiation. With some of the 
heavier molecules it was necessary to heat the 
cell. This presented a difficulty, because the 
windows which were rocksalt plates were always 
a little cooler than the cell walls around which 
the heating coil was wound, and they collected 
a film of liquid. This was avoided by putting 
the material to be studied in a little cup or 
cistern outside the cell but opening into it. This 
cistern was five or ten degrees cooler than the 
cell which could be brought up to 150°C. With 
those materials which required the heated cell, 
of which dicyclopentadiene was one, the deter- 
mination of the density of the vapor in the path 
of the radiation becomes an estimate, based on 
the boiling point of the material and on what is 
known about its change of vapor pressure with 
temperature. Unless otherwise specified on the 
charts, the absorption occurred in a cell 12 cm 
long. 

A good deal of work has been done on the 
infrared spectra of hydrocarbons.’ This work, 


6 Schaefer and Matossi, in Das Ultrarote Spektrum, in 
Article 33 (page 266) on the spectra of organic liquids and 
vapors, give a summary of the infrared absorption spectra 
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however, presents no systematic examination of 
the types of compounds present in gasoline, and 
with the exceptions noted the spectra are carried 
only through the very near infrared region. In 
fact, the work with relatively high dispersion is 
concerned largely with the absorption bands near 
3u, which seem to be shown by all hydrocarbons, 
and are taken to arise from variations in the 
single hydrogen-carbon separation. The charts 
given with this paper, while not characterized 
by high resolution, present the spectra of certain 
substances in the important region between 1 
and 30u. The charts themselves do not extend 
to 30u, but no absorption bands (at wave-length 


of the hydrocarbons as known up to the end of 1929. The 
most important work there cited on the present problem is 
that of Coblentz. J. W. Ellis also, Phys. Rev. 27, 298 (1926), 
discusses the work of Coblentz and adds much to the known 
spectra of liquid materials, including benzene, toluene and 
certain chain compounds in the region up to 2.8u, which is 
accessible with a quartz spectrograph, but he gives no new 
material for gases, nor for the longer waves. Higher 
resolution has been applied to certain substances in limited 
regions of the spectrum by Ellis, by Meyer, Bronk and 
Levin, by Brackett, and by Bell. Work published since 1929 
includes the following: 

R. B. Barnes, Infrared Absorption of Some Organic 
Liquids under High Resolution. Phys. Rev. 36, 296 (1930). 
Water, benzene and toluene are studied in the region 3 to 
4u. Also in Phys. Rev. 35, 1524 (1930), Barnes describes 
work on benzene and eight other cyclic compounds in 
liquid form in the region 3 to 3.6u. 

J. W. Ellis, Faraday Society Transactions 25, 611 (1929). 
These are reports at the symposium on molecular spectra 
and molecular structure held by the Faraday Society. 
The work of Ellis (page 888) concerns liquids and the 
spectral region between 1 and 3u. 

Included in the extensive literature of the Raman 
spectrum are recent articles of interest in connection with 
the molecules here studied. It is sufficient to cite two: 

P. Krishnamurti, Complete Raman Spectrum in Relation 
to Infrared Absorption. Part II. Benzene, Cyclohexane and 
Octane. Indian J. Phys. VI, 543 (1932). 

E. Canals, M. Godchot and G. Canquil. Raman Spectrum 
of Some Cyclic Hydrocarbons. Comptes Rendus 194, 1374 
(1932). 

Not a great deal of work bearing directly on the present 
problem has appeared from the time that the observations 
given in this paper were completed (in October, 1931) until 
the date of this writing (July, 1932). However, a paper by 
P. Lambert and J. Lecomte, Comptes Rendus 194, 960 
(March 14, 1932), gives the location of bands of cyclohexane 
and cyclohexene and certain heavier compounds. These 
bands check in wave-length with those of cyclohexane and 
cyclohexene given here. No curves are shown, however, nor 
is any indication given of the shape or nature of the bands. 


less than 30u) were found beyond the limits 
given, and if the hydrocarbons have any such 
absorption bands they are very weak. These 
spectra constitute a first step in the study of the 
molecular structure of the compounds considered, 
a study which might be carried further by the 
use of higher resolution. It is possible to compute 
in some cases what resolution will probably be 
necessary in order to observe the fine structure 
of the absorption bands. In the spectrum of 
cyclohexane, for example, there is a three part 
band with center of 14.854 which has distinctly 
the appearance of a vibration-rotation band 
with a zero branch. The separation of the peaks 
of the P and R-branches, if such they are, is 
about 0.6u so that 


Av= (3 X10" X 0.6) /15?=8 x 10" 


vibrations per second. Assuming that the mo- 
ment of inertia concerned in the production of 
this band is one of the two equal moments of 
inertia for a thin disk, it comes out from simple 
kinetic theory to be about 710-*° g cm*®. A 
better calculation as made by Dennison in work 
not yet published gives J=14X10-** g cm’. 
By using the latter value, the fine structure 
separation is expressed by the value of the 
frequency interval between successive lines as 
Av=1.2 X10", or about 0.4 waves per cm. This 
is about 1/70 of the doublet separation. The 
very similar band of benzene at 9.54, where Ad 
is 0.25, gives almost the same moment of inertia. 
This band has been worked over with a combi- 
nation prism and grating spectrometer. The 
theoretical resolving power of this spectrometer 
was such that the separation of the band into 
its combination lines appeared to be possible, 
but careful trials did not reveal any fine struc- 
ture. If we assume the cyclohexane molecule to 
be essentially flat, the value 27 for the moment 
of inertia about an axis perpendicular to the 
plane gives the distance between neighboring 
carbon atoms as 1.4 10-* cm. 

The spectra shown in Figs. 1-25 suggest the 
possibility of a beginning in qualitative analysis 
of hydrocarbon mixtures by means of infrared 
spectra. Nearly every one of the gasoline compo- 
nents has some characteristic line or lines not 
shown by any other component. The spectra of 
benzene, cyclohexane and cyclohexene, for ex- 
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ample, though similar in certain respects, show 
unmistakable differences. That of cyclohexene 
has unsymmetrical bands, the location of which 
could hardly be predicted from the wave-lengths 
shown by the absorption of benzene and cyclo- 
hexane. In this connection something may be 
said about the sources and the purity of the 
materials employed. All the samples were sent 
to the university from the Research Laboratory 
of the General Motors Corporation at Detroit. 
Mr. Paul M. Cramer of that laboratory has 
prepared the following statement about the 
materials: 


With the exception of nickel carbonyl, germanium 
tetraethyl, tetramethyl tin and methyl nitrate, all the 
compounds included in this work were either prepared 
or repurified in this laboratory. Many of the com- 
pounds were obtained from outside sources. In only 
one or two cases is there any reasonable doubt about 
the purity of the compounds. The boiling-point range 
of the compounds was within narrow limits, 0.5°C or 
less, and the physical constants were in close agreement 
with the accepted values given in the literature. 

n-Hexane, benzene, cyclohexene, cyclohexane, 
toluene, methyl cyclohexane, n-pentane, nickel 
carbonyl, methyl nitrate and methyl iodide were 
obtained from the Eastman Kodak Company. 

Diiso butylene and 2, 2, 4-trimethyl pentane were 
obtained from the Ethyl Gasoline Corporation. The 
sample of diiso butylene consisted of a mixture of two 
isomers, 2, 2, 4-trimethyl pentene-4 and 2, 2, 4- 
trimethyl pentene-3. 


Dicylcopentadiene, a coal-tar product obtained 
from the Gesellschaft fiir Teerverwertung, was 
carefully fractionated under reduced pressure. Since 
the reaction of cyclopentadiene to dicyclopentadiene is 
reversible at 170°C, there are probably traces of 
cyclopentadiene and also higher polymers at room 
temperature. 

The aniline was purified by the recrystallization of 
the hydrochloride prepared from aniline obtained from 
the Dow Chemical Company. 

The sample of tetraethyl germanium was furnished 
by L. M. Dennis of Cornell University. 

The tetramethyl tin was obtained from C. A. 
Kraus of Brown University. 


The work at the university was done with the 
assistance of Mr. P. E. Martin, Mr. John Strong 
and Mr. V. H. Fraenckel. It has profited by the 
cooperation of many members of the staff of 
the Physics Department, and in particular by 
the help of Professor Barker and Professor 
Randall. 

Some points of interest appear in connection 
with the spectra of the compounds containing 
the tetra-valent metals, tin, lead and germanium. 
They show no double or triple bands, such as 
appear in the spectra of benzene and cyclohexane, 
but they have, like the hydrocarbons, strong 
bands at about 3.4 and in the neighborhood of 
7u. These materials require the use of the heated 
cell, and these are the only spectra not carried 
beyond the range of the Wadsworth spectrom- 
eter, that is beyond 16x. 
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General Properties of Steady Temperatures in Solids Partially Exposed to Gas 


R. V. Cuurcnit., Department of Mathematics, University of Michigan 
(Received January 3, 1933) 


Five general properties of steady temperatures in solids 
of arbitrary shape are derived here. Although the deriva- 
tions are mathematical they depend almost solely on the 
very fundamentals of the mathematical theory of heat 
conduction. The solid has any portion of its surface held at 
a fixed temperature distribution and the rest exposed to 
gas. The rate of transfer of heat between gas and surface 
is assumed to be any function of the position of the point 
and the temperature of the surface, for a given state of the 
gas. For a fixed temperature this function is called the 
emissivity. Briefly stated, the important properties are: I. 
In two geometrically similar solids exposed to the same 
surface temperature conditions at corresponding points, 
except that their emissivities are inversely proportional to 
their sizes, the temperatures at corresponding points are 


equal. I(a). If the entire surfaces of two geometrically 
similar solids are given fixed temperature distributions the 
temperatures at corresponding points are equal everywhere 
if they are equal on the surfaces. II. An increase in the size 
of a solid causes the same temperature change as an 
increase in its emissivity. IV. An increase in the size of a 
solid causes an increase in the temperature at each point 
if the solid is receiving heat from the gas, and a decrease 
if the solid is losing heat to the gas. V. From Newton’s law 
of surface heat transfer there follows a linear relation 
between the steady temperature at any point and the 
uniform temperatures of the gas and the fixed-temperature 
surface. These properties are discussed at the end of the 
paper. Some of the results of IV seem contrary to expected 
results. 


N a mathematical study of the thermal design 


of spark plug insulators' the author found 
several properties of temperatures in solids 
partially exposed to gas which could be general- 
ized so as to be applicable to a variety of prob- 
lems in physics and engineering. These general- 
ized properties, which are believed to be new, 
are derived here from the fundamentals of the 
mathematical theory of heat conduction. They 
are of special importance in problems in which 
solids of the same shape but constructed to 
different scales are subjected to similar temper- 
ature conditions (see Properties IV and I(a)). 

Except in Property V, the law assumed for 
the transfer of heat between the gas and the 
solid is very general, so it seems that few re- 
strictions on the usefulness of the properties 
should be experienced on its account. 

All temperatures considered here are those 
which exist after a steady state of temperatures 
is reached. The solids are assumed to be homo- 
geneous and isotropic, and to have constant 


1 This work was done for the Department of Engineering 
Research, University of Michigan. To the AC Spark Plug 
Company of Flint, Michigan, which company suggested 
and supported this work, the author makes grateful 
acknowledgment. 
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thermal conductivities. Their shapes are arbi- 
trary. 

The properties derived here apply to a solid 
exposed to a gas over any regions of its surface. 
To the remainder of the surface any definite 
temperature distribution is assigned. These two 
portions of the boundary of the solid are referred 
to hereafter as the exposed surface, g(x, y, z)=0, 
and the fixed-temperature surface f(x, y, z)=0. 

At the exposed surface, for any given state of 
the gas, a general relation is assumed between 
the coordinates of the point, the temperature T, 
and the temperature gradient normal to the 
surface at the point: 


dT /dn= ¢(x, y,2,T) on g(x, y,2)=0, (1) 


where 7 is the outward drawn normal to the 
surface. 

In the first four properties this general law 
(1) is restricted only to this extent: It is assumed 
that the function (x, y,z, 7) cannot increase 
with T at any point. This is a natural restriction 
since it implies that the absolute value of the 
temperature gradient normal to the surface, or 
of the rate of transfer of heat across the surface, 
cannot increase when the surface temperature 
approaches the gas temperature. In Property V 
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a very special case of Eq. (1), known as Newton’s 
law of surface heat transfer, is used; and in 
Property I(a) no gas is involved. 

The rate of transfer of heat across the exposed 
surface depends upon the value of the function 
¢. For a fixed T we shall refer to this function 
as the emissivity function. Then by an increase 
in this function, or an increase in emissivity, 
we shall mean a change to a new function 9’ 
which is greater than @ at each point of the 
exposed surface, for the same 7; that is, 
¢'(x, ¥, 2, T) > (x, y, 2, T) for all T, x, y, 2. 


PROPERTY I 


Consider any two geometrically similar solids 
with corresponding portions of their boundary 
surfaces held at fixed temperatures and the re- 
maining portions exposed to gas. If the tempera- 
tures are the same at corresponding points of the 
fixed-temperature surfaces and if the emissivities 
at corresponding points of the exposed surfaces are 
inversely proportional to the sizes of the solids, 
then the temperatures at corresponding points 
throughout the two solids are equal. 

To prove this property we must first make 
clear the meaning of similarity. Two solids are 
geometrically similar if their boundary surfaces 
can be placed relative to each other so that they 
are in perspective with respect to some point O. 
That is, if P is any point on one surface then 
there is a point P’ on the other so that the 
vectors (OP’) and (OP) satisfy the condition 
(OP’)=a(OP), where a@ is the same for all 
points P. 

Then if O is taken as the origin a point 
(x, y, 2) on the first surface has (ax, ay, az) for 
its corresponding point on the second; distances 
between corresponding pairs of points have the 
ratio 1/a, while the directions of the lines joining 
these pairs are the same. It also follows that if 
S(x, y,z)=0 is the equation of the boundary 
surface of the first solid then S(x/a, y/a, z/a)=0 
is the equation of that of the second, and the 
normals to these two surfaces at corresponding 
points are parallel. 

The differential equation and boundary con- 
ditions which determine the temperature in the 
first solid are 


#T /dx°+ /dy?+ eT /d2=0 


throughout the solid, (2) 


T=T)(x, y,z) on the surface f(x, y,z)=0, (3) 


dT /dn= $(x, y, 2, T) 


on the surface g(x, y, z)=0, 


(1) 


where 7;(x, y, z) is any function of x, y and z. 
Let T= F(x, y, 2) be the temperature distribution 
determined by these conditions. 

Let us substitute the new independent vari- 
ables x’=ax, y'=ay, 2’=az throughout this 
problem, where a is a constant. The solution 
becomes T= F(x’ /a, y’/a, 2’/a). If u, v are the 
direction cosines of the normal to the surface 
g(x, y, z)=0 at (x, y, 2), then 


dT oT oT oT oT 
dn oy dz ax’ ay’ 


Since the direction cosines of the normal to 
g(x, y, z)=0 at (x, y, 2) are the same as those of 
the normal to g(x’/a, y’/a, 2’/a)=0 at (x’, y’, 2’) 


the expression in parentheses in the last equation 
becomes dT/dn’, when n’ is the normal to 
g(x’ /a, y’/a, 2’/a)=0. 

In terms of the new variables then, Eqs. (2), 
(3), and (1) become 


#T 
for all x’, y’, 2’, (2°) 
T= T(x’ /a, y’/a, 2’ /a) 
when f(x’/a, y’/a, 2’/a)=0, 
dT /dn'= /a, y’/a, 2’/a, T) 


when g(x’/a, y’/a, 2’/a)=0. 


(3’) 


(1) 

When we interpret x’, y’, z’ as the coordinates 
of a point, the surfaces f(x’/a, y’/a, s’/a)=0 and 
g(x’/a, y’/a, 2’/a)=0 are in perspective with the 
surfaces f(x, y,z)=0 and g(x, y,s)=0. Hence 
Eqs. (1’), (2’), (3’) are the complete set of 
conditions for the problem of temperature 
distribution in a solid geometrically similar to 
the first and a@ times as large. 

According to Eq. (3’) the temperature at 
(x’, y’, 2’) on the fixed-temperature surface is 
the same as that at the corresponding point of 
the first solid. Eq. (1’) shows that for equal 
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temperatures at corresponding points of the 
exposed surfaces of the two solids the normal 
gradient of 7 at the exposed surface of the 
second is 1/a times its value at the corresponding 
point of the first; that is, the emissivity of the 
second is 1/a times that of the first. 

Since the function T= F(x’ /a, y’/a, 2’/a) satis- 
fies the conditions (1’), (2’), (3’), it gives the 
temperature distribution in the second solid. It 
shows that the temperature at (x’, y’,2’) is 
the same as that at the corresponding point 
(x’/a, y’/a, 2’/a) in the first solid. Hence Prop- 
erty I is proved. 


Ppeorerty II 


An increase in the size of a solid causes the 
same temperature change as a proportionate in- 
crease in its emissivity. 

More precisely, consider a given solid and a 
second one geometrically similar to it but a 
times as large. Let these solids have equal 
temperatures at corresponding points of their 
fixed-temperature surfaces, and equal values for 
their emissivity functions at corresponding points 
of their exposed surfaces, that is, $’(x’, y’, 2’, T) 
= o(x, y, 2, 7). Then the temperature at each 
point of a third solid, identical to the first except 
that its emissivity function is @ times as large, 
will be the same as the temperature at the 
corresponding point of the second. 

This is a corollary of Property I. Let 
(x, y, z, T) be the emissivity function at (x, y, 2) 
of the first solid. Then at (x’, y’, 2’) of the second, 
where x’=ax, y’=ay, 2/=az, the emissivity 
function is $(x’/a, y’/a, 2’/a, T), and at (x, y, 2) 
of the third it is a¢(x, y, z, T). The emissivities 
of the second and third are therefore inversely 
proportional to their sizes so the temperatures 
at corresponding points are equal, according to 
Property I. 

In the next two properties the assumption 
will be made that the temperature determined 
by our conditions (1), (2), (3) is continuous over 
the surface of the solid as well as throughout 
its interior. 


Property III 


Except at points on the fixed-tem perature surface, 
the temperature at each point of a solid increases 
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as a result of an increase in the value of its emis- 
sivity function at all points of the exposed surface, 

Let us state this property for each of the two 
cases ¢>0 and ¢<0. In a solid which receives 
heat from a gas at all points of its exposed 
surface the temperature at any point cannot 
decrease when the value of the emissivity 
function is increased. But in a solid which loses 
heat to a gas at all points of its exposed surface 
the temperature at any point cannot increase 
when the absolute value of the emissivity 
function increases. It is only on the fixed- 
temperature surface that the temperature re- 
mains unchanged. 

This property may seem too obvious physically 
to deserve a mathematical proof; but it leads 
directly to Property IV which, in some of its 
applications, seems to defy intuition. 

In a solid with emissivity function @ the 
temperature T satisfies Eqs. (1), (2), (3). Let T’ 
be the temperature in a solid of the same size 
and shape subject to the same surface tempera- 
ture conditions except that its emissivity function 
is increased : 


¢'(x, y, 2, T’) > y, 2, T’) on g(x, y, 2) =0. (4) 


This increase may be due, for example, to a 
change in the state of the gas. 
Then 7” satisfies the conditions 


throughout the solid, 
T’=T,(x, y,z) on f(x, y, z)=0, 


dT’ /dn= 9'(x, y,2, 7’) on g(x, y, 2)=0. 


Now let u=7’—T; then according to the 
above conditions on 7 and 7’, u must satisfy 
the equations 


/dx? + 0u/dy?+ Pu/dz?=0 
throughout the solid, (5) 


u=0 on f(x, y, 2)=0, (6) 


du/dn= ¢'(x, y, 2, T’) — o(x, y, 2, T) 
on g(x, y,2z)=0. (7) 


Suppose that the function « which satisfies ' 
these three conditions has at some point P a 
minimum value which is negative. Then the 
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point P must lie on the surface g(x, y, z)=0. 
For u is a harmonic function, according to Eq. 
(5), and the minimum value of such a function 
is assumed only on the boundary of the solid,’ 
unless u is the same at all points. However, u 
is not negative on f(x, y, z)=0 because of Eq. (6). 

Hence at some point P of g(x, y,z)=0 we 
should have 7 >T’ and, since ¢ does not increase 
with 7, it would follow that 


y, 2, (x, y, 2, T). 
Condition (4) then gives the inequality 
(x, y, 2, T’) > o(x, y, 2, T’ )=o(x, y, 2, T) at P, 


and according to Eq. (7), this means that du/dn 
is positive at P. But this implies that wu increases 
toward its minimum value at P, which is 
obviously a contradiction. The minimum value 
of u is therefore zero. 

Now u cannot be the same at all points since 
this would mean u=0, and T=T”’, everywhere, 
whereas the conditions (4) and (7) require du/dn 
to be different from zero on g(x, y, z)=0. More- 
over “% cannot assume this minimum value on 
the exposed surface since the conditions (4) and 
(7) require du/dn to be positive there. Hence 
u=( on the fixed-temperature surface only, and 
u>0 elsewhere. 

The temperature 7” is therefore greater than 
T at each point not on the fixed-temperature 
surface, where 7’= 7, and Property III is proved. 


Property IV 


An increase in the size of a solid to one geo- 
metrically similar to the original causes an increase 
in its temperature if the solid is receiving heat from 
the gas, and a decrease in its temperature if the 
solid is losing heat to the gas. 

To make this important property clear, con- 
sider two solids which are geometrically similar 
and subject to the same fixed temperatures and 
equal emissivity functions at corresponding 
points of their surfaces; that is, J (x, y, 2) 
= T)'(ax, ay, az) on fixed temperature surfaces, 
o(x, y, z, T)= ¢’(ax, ay, az, T) on exposed sur- 


faces. Then the temperature at each point of the 


* See Riemann-Weber’s Differentialgleichungen der Physik 
I, 1925, p. 571, or MacMillan, Theory of the Potential, 1930, 
p. 134. 
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larger solid is higher than that at the corre- 
sponding point of the smaller provided heat is 
flowing from the gas to the solids. If heat is 
flowing from the solids to the gas then the 
temperature at each point of the larger solid is 
lower than that at the corresponding point of 
the smaller. At corresponding points of the 
fixed-temperature surfaces, of course, the temper- 
atures are equal. 

This property follows easily from Properties 
II and III. A solid which is @ times as large as 
the original but subject to the same surface 
temperature conditions has a temperature at 
each point equal to that which would occur at 
the corresponding point of the original if the 
emissivity of the original were multiplied by a 
(Property II). Take a>1, and suppose that 
heat flows from gas to solids (¢>0 everywhere). 
It follows from Property III that the tempera- 
ture at each point of the large solid, except for 
the fixed-temperature surface, is higher than 
that at the corresponding point of the original. 
Similarly if we suppose that heat flows from the 
solids to the gas (¢<0 everywhere), it follows 
from Property III that the temperatures in the 
large solid are lower than the corresponding 
temperatures in the original. 


PROPERTY V 


When Newton's law of heat transfer at the 
exposed surface applies and the temperature Ty is 
uniform over the fixed-temperature surface, the 
steady temperature T at any point of the solid 
bears a linear relation to the temperature of the gas 
T, and the temperature Ty: 


T= WT,+(1—W)T. (8) 


The quantity W depends for its value upon 
the design of the solid, the coordinates of the 
point at which 7 is measured, and the thermal 
coefficient in Newton’s law (10) below. 


The temperature 7 is determined by the three 
conditions 


/dx? + 0°T /dy?+ /d2’=0 at all points, 
T=T,) on f(x, y,2)=0, 
dT/dn=h(T,—T) on g(x,y,2)=0, (9) 


where h is a constant. If we substitute the new 
variable 
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W= (T—To)/(T,— To) (10) 


these three conditions become 0? W/dx?+ 0°?W/dy? 
W=0 on f(x, y,2)=0, dW/dn 
=h(1—W) on g(x, y, z)=0. 

But these are the complete set of conditions 
for determining the steady temperature W at 
any point of the given solid whose surface 
f(x,y, z)=0 is held at zero temperature while 
the remainder of its boundary g(x, y,z)=0 is 
exposed to Newton’s transfer of heat from a gas 
with unit temperature. Hence W depends upon 
x,y,2, the value of the constant h, and the 
design of the solid. Since the linear Eq. (9) is 
another form of Eq. (11), Property V is proved. 


DISCUSSION OF PROPERTIES 


It should be noted that insulated regions of 
the boundary surface are not excluded in the 
above properties since ¢=0 on the region 
insulated is the mathematical condition for this 
case. In such cases a minor change is necessary 
in Property IV, however. In our proof, extended 
by assuming a continuously varying T as ¢ 
approaches zero over any region, we would not 
exclude the possibility of T being the same at 
corresponding points on these insulated surfaces 
of the two solids. 

When the two solids of Property I are placed 
in perspective it is clear that the isothermal 
surfaces in the two solids for any given tempera- 
ture are also in perspective. This is true for any 
two geometrically similar solids with their entire 
boundary surfaces held at corresponding fixed 
temperature distributions. To prove this all we 
need do is to modify the conditions in our proof 
of Property I as follows. Let f(x, y, z)=0 be the 
entire boundary surface of the first solid,: so 
that f(x’/a, y’/a, 2’/«)=0 is that of the second, 
and drop the conditions (1) and (1’), so that 
the function ¢ is not involved. The following 
important special case of Property I is thus 
obtained: 


I(a) 


If two geometrically similar solids have their 
entire boundaries held at fixed temperature distribu- 
tions such that temperatures at corresponding points 
of the two boundaries are equal then the tempera- 


tures at corresponding points throughout the two 
solids are equal. When the two solids are placed in 
perspective their two isothermal surfaces corre- 
sponding to any given temperature are in pers pec- 
tive, with the same ratio a of perspectivity as the 
boundaries. 

This raises the interesting question of the 
relation between corresponding isothermal sur- 
faces in geometrically similar solids subject to 
more general surface temperature conditions. 
It is just recently that the examination of the 
nature of families of isotherms for variable 
temperatures was begun.’ 

It is important to note that the proof of 
Property IV is independent of the shape of the 
solid and the amount of its surface area which is 
exposed to the gas. Suppose, for example, that 
a slender solid cylinder with length L and 
diameter D has its entire surface, except for one 
base, held at a fixed high temperature, while the 
other base is exposed to cool air at fixed temper- 
ature. Consider a second cylinder with dimen- 
sions L’= 5L, D’=5D subject to the same surface 
temperature conditions. Let P; be a point in the 
small cylinder at a small distance c from the 
lateral surface and at an equal distance from 
that base whose temperature is fixed. In the 
large cylinder the corresponding point Pz» is at a 


- distance 5c from the corresponding surfaces. 


According to Property IV, the steady tempera- 
ture at P2 is lower than that at P;. 

This result becomes even less evident intui- 
tively if the cylinders are replaced by two 
similar truncated cones with only their small 
bases exposed to cool air. Other more complex 
shapes illustrate still better the possible devi- 
ations of the results of Property IV from the 
expected results. 

In Property V, h=E/K where E is the 
external conductivity and K is the thermal 
conductivity. Hence if Newton’s law (10) 
applies, Properties I, II, III can be stated 
easily in terms of K for a fixed E as well as in 
terms of the emissivity h. 

The linear relation (9) is a convenient one to 


3E. Kasner, Geometry of the Heat Equation, Proc. Natl. 
Acad. Sci. 18, 475-480 (1932). Also see abstracts of papers 
by Kasner (p. 811), Comenetz (p. 803), and Fialkow (p. 
803), in the Bulletin of the American Mathematical 
Society, 38 (1932). 
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use to determine the applicability of Newton’s 
law in a physical problem. If measurements of 
T at a fixed point are plotted against measured 
values of 7, with Jo fixed, or against values of 
T, with T, fixed, the deviation of these points 
from a straight line means deviation from 
Newton’s law. 

It is evident that two sets of corresponding 
values of 7, and T, where T is measured at a 
fixed point, say the hottest point, determine the 
values of ZT) and W. Thus in a problem of 
design the values of Ty) and W may be determined 


by the extreme values allowable for T and T,,. 
A body which has this value for W at its hottest 
point may be obtained by adjusting either the 
size, shape, or thermal coefficients. According 
to Properties II and IV however, any value of 
the relative temperature W made possible by 
adjusting thermal coefficients can also be ob- 
tained by adjusting the size of the body. That 
is, for given extreme allowable temperatures and 
given thermal coefficients there is in general an 
optimum size for any shape rather than an 
optimum shape. 
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Stress Distribution in Wedges with Arbitrary Boundary Forces* 


J. H. A. Brantz, Aeronautics Department, California Institute of Technology 
(Received December 20, 1932) 


Stress distribution under plain strain in a corner of any 
angular magnitude. A generalization of the method of von 
Karman, applicable to the wedge. for any force distribution 
over the straight boundaries is outlined. In the first 
application the stress function is found for a concentrated 
force acting at any point of the boundaries of the three- 
quarter plane. Certain stresses are plotted. The stress 
functions are given without mathematical details for the 
following boundary conditions in the three-quarter plane: 


(a) Uniform tension over a part of the boundary; (b) 
linear distribution over a part of the boundary; (c) 
superposition of (a) and (b), giving hydrostatic distri- 
bution, with a plot of certain stresses. The discussion 
contains the stress functions for uniform continuous tension 
on one or both boundaries and points out the very inter- 
esting paradox that stresses may be finite for certain 
continuous loadings, but become infinite if a portion of the 
load is removed. 


I. INTRODUCTION 


HE knowledge of two-dimensional stress 

distributions is useful for the solution of 
many problems in mechanical and civil engineer- 
ing, when the usual methods of strength of 
materials are not applicable or not sufficiently 
accurate. The stress distribution in the half-plane 
under the influence of a concentrated force acting 
on the straight boundary is of fundamental im- 
portance, and the solution of this problem is due 
to Boussinesq. Several extensions of his results 
have been made by Michell,' especially by means 
of the so-called method of inversion, the only con- 
formal transformation of bi-harmonic solutions 
which leaves the boundary stresses undisturbed. 
By superposition of Boussinesq’s solutions the 
stress distribution in the half-plane is easily 
obtained for an arbitrary load distribution along 
the straight boundary. 

The stress distribution in a plane strip has been 
discussed by several authors and special reference 
is made to a paper by Dr. Th. von Karman’ con- 
taining a general method for the determination 
of the stresses due to an arbitrary load along the 
boundaries. The paper is not yet available in the 
English language, therefore a brief description 
shall be given of v. Karman’s formulation of the 
problem. 


* Excerpt from Doctor s thesis by the author. 
1A, E. H. Love, Elasticity, page 216. 
2y. Karman, Aachen Abhandlungen, 1927, p. 7. 


A two-dimentional elastic system? is in equi- 
librium if the stresses o,, ¢, and 7 are the second 
partial derivatives of an arbitrary function 


F(xy): 
o,=0°F/dy"; o,=0F/ 0x"; —0°F/dxdy. (1) 


In order to be compatible with the stress-strain 
relations in accordance with Hooke’s law, F must 
be a solution of the equation 


O'F/dx*+ 20°F /dx*dy?+ 
or in Gibbs’ notation 
ViF=0. (2) 


F, known as Airy’s function, and its first 
partial derivatives can be interpreted in terms 
of the quantities usually employed and defined 
in “strength of materials,’’ namely, shear and 
moment of the external forces. By integrating 
(1) along the boundary between points A and B 
it is seen that the differences 


and 


represent respectively the resultant Y and X 
components of the load on the boundary between 
A and B and, by integrating once more, it is seen 
that the difference (F),—(F), represents the 
increase in the bending moment, in going from 
A to 


3 See Apendix I. 
‘ See: Riemann-Weber, Differentialgleichungen der Physik, 
Il. 
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Hence, if all loads are assumed applied at the 
upper surface of the strip, x-axis (see Fig. 1), the 
usual moment curve becomes identical with 
F(x, 0) for the strip in question and the problem 
of stress distribution becomes simply the deter- 
mination of F(x, y) for given boundary values. 


Y 
A 
y=O 
h 
y=-h 
Fic. 1. 


It can easily be verified that the following 
functions are particular integrals of (2): 


A,cos mx Cosh my, A;sy cos mx Cosh my, 
Azcos mx Sinh my, Ag¢y cos mx Sinh my, 
A;sin mx Cosh my, Azy sin mx Cosh my, 


A,;sin mx Sinh my, sin mx Sinh my, 


and if we consider the coefficients A, etc., func- 
tions of the parameter m the Airy’s stress func- 
tion can be written: 


f [(A, Cosh my+-A xy Cosh sey 
0 


+Agy Sinh my) cos mx+(A; Cosh my 
+A,Sinh my+Azy Cosh my 
+Asy Sinh my) sin mx]dm, (3) 


where Aj, etc., are to be determined by the 
boundary conditions. 

If the depth of the beam is h, the boundary 
conditions to be satisfied by F are: dF/dy=0, 
F= M(x) at the upper boundary, y= 0;0F/dy=0, 
F=(0 at the lower boundary, y= —A. As in Fig. 1, 
the moment M(x) is assumed to be zero except 
in the region x= a, to x=dz2 and we now express 
M(x) as a Fourier integral: 


M,= (a cos sin mx)dm, 
0 


where a and } are determined by 


1 
=-]| M(a)-cos mada, 


1 
AM (a) sin mada. 

Te a, 
We now form the four boundary equations which 
must hold for all values of x. Hence, by equating 
to zero the coefficients of cos mx and sin mx we 
obtain eight other equations which determine the 
coefficients A, etc. (later in this paper this 
process is carried out in detail), and by the sub- 
stitution of these in (3) we finally obtain the 
stress function F. By (1) the stresses are then 
found. 

In other words, if it is possible to represent 
the moment distribution M(x) as a Fourier’s 
integral the solution of the stress problem is 
immediately obtained. 

v. Karman and F. Seewald have applied this 
method to several important cases and obtained 
interesting results, among others the “‘corrected”’ 
relation between the bending moment and shear 
at a point of a beam and the curvature of the 
central line 


(A) 


where K denotes the curvature, E and G Young’s 
modulus and the corresponding modulus in 
shear, Poisson’s ratio, A area of cross section, 
I moment of inertia, M and P the bending 
moment and shear at the point in question. The 
second term corresponds to the influence of 
shear on the deflection. 

A new and interesting conclusion of this theory 
is that the deflection of the central line of a long 
beam is wave shaped outside the loaded region. 


II. SOLUTIONS FOR ARBITRARY BOUNDARY LOADS 
ON A WEDGE 


The method described above can be applied to 
any angular portion (vy) of the plane with any 
distribution of both normal and tangential forces 
over the straight boundaries. 

Introducing polar coordinates and stress func- 
tion must satisfy the equation 


VtF(r, 0) = (5) 

with corresponding stresses, see Fig. 2: 
(6) 
op= 0 F/dr’, (7) 
t= —(0/dr)(0F/rd8). (8) 
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For the present purpose we will employ the 
functions obtained by taking real and imaginary 


parts of 
(x and y)-(2'™+2-™), (9) 


where m is real and to be determined by the 
boundary conditions. 


We obtain the following particular solutions 
of V'F=0: 


¢;=xcosmétCosh m8, 
¢g2=xcosméSinh mé, 


¢s= ycosméCosh 
¢s= ycosméSinh 
ysinméSinh mé, 
¢s= ysin mé Cosh m8, 


(9a) 


¢g3=xsinméSinhmé, 
¢4=xsinméCosh 


where §=log r. 

The coefficients A, are now defined as func- 
tions of the parameter m. The stress function 
may then be expressed as 


F=> 


n=1e/9 


A n(™) dm 


0 


from which the stresses are obtained by (6), etc. 

Two typical cases of boundary loads must be 
considered : 

(a) The load may be continuous (not necessarily 
uniform) over portions of the straight boundaries. 
Let it be assumed that the load is distributed 
over finite portions of the boundaries or, if the 
load extends to infinity, that the load per unit 
length becomes zero at infinity. Let the bound- 
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aries be the lines 02=0 and 6@=y, and let the 
normal loads be respectively po and p, and the 
tangential loads be go and q,. These quantities 
are functions of 7, measured from the inter- 
section of the boundaries, or functions of = log r. 
Using Fourier’s double integrals, we can write 
in the following form®: 


pot) = [cos me f mé "pala cos mada 


+sin mE] po(a)-sin (11) 


where the load po(£) is assumed to exist only 
between the points £; and £. Similar expressions 
are obtained for p,, go and gy. Therefore, the 
boundary conditions are 


Po(t) expressed asin (11) (12) 
(13) 


(15) 


(7) Yo, 


(6) Pr, (14) ys 


where go» and 7 are given by (7) and (8). 

By now equating to zero the coefficients of 
cos mé and sin mé in the above equations, eight 
other equations are obtained which determine 
the eight functions A(m) and consequently the 
stress function F by substituting A,(m) in (10). 

(b) In the case of concentrated loads the 
boundary conditions are obtained in the fol- 
lowing manner: It was already shown that 
the differences and 
represent, respectively, the 
Y and X components of all the boundary forces 
between A and B. Therefore, in the case of con- 
centrated forces, these differences become dis- 
continuous (similar to the shear polygon for 
concentrated forces) and the jump in these 
values at a certain point is equal to the re- 
spective components of the concentrated force 
at that point. Consider for instance unit normal 
force acting at point r=a of the boundary @=0. 
Obviously where « 
is a small positive quantity. If there is no con- 
tinuously distributed load on the boundary 
6=0, then by (12) (0? F/dr*)o= 

Hence dF dr is a constant on either side of 


5 See Fourier Sertes by W. E. Byerly. 
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r=a. These conditions can all be satisfied by 


writing® 
OF 1 i 
0 


Or m 


This integral has the value — } for &=log (r/a) <0 
or r<a and +} for &>0 or r>a. The integral’: 


1 sin mE 

am= |} when (17) 
0 m? +n? Owhen 

with its ramifications is often very useful in 

expressing boundary functions on Fourier in- 

tegral form. (See Fig. 3.) 


III. ApPLIcATION 


In the first application of the foregoing, the 
stress function F corresponding to a unit normal 
traction at point r=a of boundary @=0 will be 
found for a reentrant right corner, i.e., y= 37/2. 

The boundary conditions are when 6=0: 
dF/dr=0 when ¢=log (r/a) <0, 0F/dr=1 when 
dF/a0=0; when F=0, aF/a0=0. 
Let F= F;. Fo is determined by the boundary 


“conditions 


(0F,/dr)>= —} for = +3 for E>0 
(OF, 0, (Fo),=9, (0Fy/00),=0. 
Hence by Eq. (16), 


In order to remove the concentrated force in- 
troduced at the origin (because of the discontin- 
uity in the first derivative) F, is determined so 
that (0F,/dr)o= for — 
(F:),= (0F,/00),=0. F; is easily found by em- 
ploying the simple solution to V4F,=0: 


F\= cos 0+ Bér sin 6 
+Crcos0+Drsin 6. (18) 


Hence, by the boundary conditions for F;: 
(0F,/dr)>= C=}; (0F\/00),>= Ar+Dr=0; (Fi), 
= —(32/2)Br—Dr=0; (OF, /00),=(32%/2)Ar 


*See Peirce, Short Tables of Integrals, No. 484. 
7 Riemann-Weber, Differentialgleichungen der Physik, I, 
page 157. 
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— Br+Cr=0. Solving these and nests in 
(18) we obtain 


F\= cos sin 6). (19) 


The first degree terms have been omitted as they 
contribute no stress. 

In order to determine Fy we employ Eq. (10) 
and formulate the given boundary conditions 
by differentiating the functions given by (9a) 
and obtain: 


OF» 

=- —dm, 
m 

(— = Ass 0, 

rd0 

—AsCs 

OF» 

rae 

]dm=0, 


where C=Cosh ym, S=Sinh ym, c=cos mé, 
s=sin mé. By equating the coefficients of c and 
s to zero we obtain the eight equations: 


A,+mA,=0; 
mAe+A;=0; mA3;+As=0; CA;+SAs=0; 
— SA;—CAs=0; CA,+ SA2—mSA;—mCA,=0; 

and solving: 

A,= —m(S?—m’)A; A;=—m’CSA; 

+CSmA; 

A;= —CSA; A;= —mC°A; 

A,y=+(S—m’)A;  As=+mCSA; 


where A= 1/2m(1+m?)(Sinh? ym —m?). 

The stress function is then found simply by 
substituting the A’s back into Eq. (10). The 
stress-function in the ?-plane for unit tensile 
force concentrated at the point (r, 6)=(a, 0) 
becomes after some reduction F= 
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F= Cosh mé—(xS-+myC) Sinh m(y— cos 0+2sin 8). (20) 


The integral F may be evaluated by the method of contour integration in the complex plane,’ 
by making use of Cauchy’s theorem of residues as outlined in the appendix. 
The stress functions evaluated in this manner valid for small and large values of |£| become: 


Fot+ Fi= 


when r/a<1, 


= + 2.629(r/a)* cos 6+2 sin 6) —a cos? 


where 


(28) 

when r/a>1, (29) 

g(@)=sin (0.5450) cos @—0.545 cos (0.5450)-sin 6—0.839 sin (0.5458) sin 8, (30) 
h(@)=sin (0.9098) cos 6—0.909 cos (0.909 6) -sin 6+0.416 sin (0.9098) sin @. (31) 


It is interesting to notice that the terms in (28) 
including those of higher order are the bihar- 
monic functions of the wedge (solutions of 
V!F=0) which correspond to boundaries free of 
forces. This can also be seen by writing 


F=Ar" cos n0+ Br" cos (n—2)0 


+(Cr" sin n6+Dr" sin (n—2)0, (32) 
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where ” may be complex and finding the values 
of n which give possible solutions in case F and 
0F/00 are both to be zero along both boundaries. 
The determinant of the four boundary equations 
gives exactly the roots corresponding to the 
exponents in (28). It is, however, only possible 
to determine the ratios B/A, etc., in this way. 
The absolute values of the constants naturally 
depend upon the location and magnitude of the 
concentrated forces, and these can only be 
specified by the use of Fourier integrals as above. 

The writer is informed that Dr. H. M. 
Westergaard, Professor of Theoretical and Ap- 
plied Mechanics, University of Illinois, in a 
communication to Dr. Th. von Karman, has 
pointed out the existence of these solutions. 

Check on boundary conditions: By (28) and 
(29) it is easily verified that 


F=0 6=0 
aF/a0=Oland 
are satisfied. 

In order to investigate if all boundary forces 
are in equilibrium we will determine the stresses 
on a circle of large radius r= R (see Fig. 4). By 


differentiating (29) and only taking such terms 
that contribute: as r= 


(127 sin cos 0)/(97?—4)R, 
sin 20/R’. 


8 See for explanation Mac Robert, Functions of a Complex 


Variable, 
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Projecting on x-axis and letting Ro =: 
Y 
x= rd0 cos 0=0. 
0 
Projecting on Y-axis: 
rd0 o, sin 0=1. 


Moments about the origin: 


cos 20] = +a. 
0 


Cosh ym-Sinh rm 


It will be seen that all equilibrium conditions 
are satisfied, the only other force being unit 
tensile force at a. The stresses for large and 
small radii are now obtained from (28) and (29) 
by application of Eqs. (6), (7) and (8). It will be 
noticed that stresses become infinite at the 
origin. 


Evaluation of stresses 


If it is only desired to know the stresses at 
certain points or the distribution on a line it is 
often simpler to specialize F in (20) before 
evaluating the integral. 

For example, to obtain the distribution of 
normal stress on the Y-axis, place @= 7/2 in (20): 


r 
Pu —— 
my (1+m?)(Sinh? ym —m?) 


1 


(m cos mE—sin mé)dm-+ Fi, (33) 


m Cosh ym-Sinh xm 


or? xr 0 


Notice that F; contributes nothing. 

By (28) and (29) we obtain the correct 
stresses for small and large radii which usually 
are the most interesting, so that only a few 
additional terms of higher powers in r would be 
needed to give the stresses in the region r=a. 
These can be obtained by the method of residues 
applied to (34) as explained above. However, if 
only a plot of stresses is desired, it is simpler to 
evaluate (34) for values of 7 in the neighborhood 
of a, by more primitive methods. The simplest 
of these would perhaps be to plot the integrand 
for various values of m and r and apply Simp- 
son’s rule. This need only be done for m=1 
because for larger m, (34) can be further sim- 
plified by placing Cosh my=Sinh my= e*™/? and 
ignoring m? in the denominator. The integral then 
becomes 


1 


cos médm. 


The stresses og obtained by evaluating (34), 
are plotted in Fig. 5. The corresponding stresses 
for the half-plane, y=, are also shown. It is 
interesting to notice the effect of the additional 
quarter plane by comparing the two stress dis- 
tributions, I and II. Boussinesq’s solution for a 


Sinh? ym —m? 


cos médm. (34) 


unit traction at x=a on the x-axis gives 


(35) 


along the Y-axis, which is plotted as the curve II 
of Fig. 5. 

If several concentrated forces are applied on 
the boundary @=0, the stress functions and 
stresses are obtained by superposing the func- 
tions computed for the various values of a. 
For a concentrated force acting on the boundary 
6=37/2, the stresses are obtained from the 
above formulae by substituting (37/2) for 
6, maintaining the same coordinate system. 

It is of interest to note that the function g(@) 
given in (30) is symmetric and h(@) given in (31) 
is antisymmetric about the bisector of the angle 
y. This holds also for the corresponding functions 
for higher powers in r/a. This would mean that 
concentrated forces of equal magnitude located 
on each boundary at equal distance from the 
vertex would have stress functions (valid near 
the vertex): 


F=0.202a(r/a)'**- (28a) 


if the unit forces are both tension and 


F=0.250a(r/a)*™ - h(6) (29a) 
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Normal stresses % on Y-axis due to 
concentrated traction at r=a,@=0. 


Curve I in the - plane, 7° =F . 
Curve Zin the -plane,y =x . 


Fic. 5. 


with traction on the boundary @=0 and pressure IV. ApbITIONAL RESULTS 
on @=32/2. This shows that the stress at the 
vertex is always infinite for partial loading of the 
boundaries. 


Other types of boundary conditions for a right 
reentrant corner have been studied the results of 
which shall be stated here omitting the mathe- 
matical details. 


(a) With uniform tension p on the boundary 6=0 between points r=0 and r=a 


r@ 
4—97? 


F= ap| Cosh m6+(xS+myC) Sinh m(y—6@) ]sAdm+ [3a cos sin (36) 
0 


(b) With linear tension pr on the boundary 6=0 between the points r=0 and r=a: 


F= op| Cosh mé—(xS+myC) Sinh m(y—6) }- [me—(m?+2)s jAdm 


[32 cos 0+2 sin 6] (37) 


(c) By superposition of a and b we obtain the stress function for hydrostatic loading (tension) 


F= op| Cosh mé+(xS+myC) Sinh m(y—6)]-(mc+2s)A-dm 


r@ 
24-97? 


(3m cos sin a). (38) 
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In Fig. 6 the normal stresses on the line @= 7/2 
because of hydrostatic loading are plotted. 
Curve II shows the same stresses for similar 
loading on the half plane (y= 7). 

It will be noticed that Eq. (23) in the Appendix 
has no real roots, except m=0, when y is less 
than z. This means that in general a partial 
loading of the boundaries of a corner will 
produce finite stresses at the origin (vertex) if 
y<za, but infinite stresses if y>za, and the 
greater y is the sharper the infinity. The half- 
plane, y=7, thus is the critical case and the 
stresses at the origin are zero if it is chosen 
outside the loaded region. 

Special cases of stress distribution in an infinite 
wedge have been treated by a number of inves- 
tigators. We find J. H. Michell® who generalizes 
Boussinesq’s solution-for a concentrated force 
on the half-plane to a concentrated force at the 
vertex of an infinite wedge, M. Levy® who finds 
a solution for a linear load distribution over the 
boundaries, S. D. Carothers” who gives solutions 
for a couple at the vertex and uniform loads on 
the boundaries. 

The stress function for continuous uniform 


° A. E. H. Love, Elasticity, page 212. 


10 Proc. Roy. Soc. (Edin.) 53, 292 (1913). 


pressure on the boundaries is easily determined 


by Eq. (32). 
= —1pr* cos? 8 (39) 


corresponds to uniform pressure on the boundary 
6=0. It will be seen that the stresses are finite 
at the vertex for any value of y. We have seen 
that a partial loading produces infinite stress at 
the vertex if y>7, so that we have the curious 
result that in this case the removal of a portion 
of the load raises the stress to infinity. 

If uniform pressure acts on the boundary 
6= 37/2 the corresponding stress function is 


= —}pr’ sin? 0. (40) 


In Fig. 6 is plotted the normal stresses on the 
base, = 7/2, due to linear pressure distribution 
on the boundary, @=0. It is important to notice 
that near the vertex the base is under high 
tension. 

Levy’s solution for linear normal pressure on 
the boundary of an infinite wedge (or triangle) 
shows linear distribution of all stresses along 
any straight line through the wedge. In practice 
this result is conveniently applied to masonry 
dams. However, the application is dubious 
because the height of a dam is of the same 
order of magnitude as the base width, so that 
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great deviations from the linear stress distribu- 
tion can be expected (at any rate near the base) 
because of the influence of the unyielding, or 
rather, infinite foundation. 

Engineers are well aware of this fact, and the 
present study grew out of an attempt to deter- 
mine the stress distribution in triangular dams 
on infinite, rigid and elastic foundations by the 
Ritz’ method. However, it was found that the 
stresses in the up- and downstream faces at the 
base did not converge when more parameters 
were included in the solution. At points only 
slightly distant from the up- and downstream 
corners Ritz’ method gave reasonable results. 

The present method enables us to compute the 
stresses in dams of triangular and other shapes 
for any load and far more accurate solutions may 
be reached in this manner. 

Photoelastic experiments are at present being 
carried out at the California Institute by the 
author for the purpose of verifying the theoretical 
results of this study. It is of special interest to 
find out how closely the equation V‘F=0 holds 
for plates of finite thickness. The stress distribu- 
tion is also being investigated in plates of arbi- 
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m-f(m)+cos 1 
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trary shapes, especially applicable to masonry 0 
dams, also the stress concentrations around open- 

ings corresponding to the inspection and drainage 
galleries which usually exist in high dams and 7 
finally the stresses caused by temperature changes a 
and chemical shrinkage. 

The author wishes to express his appreciation 
to Dr. Th. von Karman, Director of the Guggen- 
heim Graduate School of Aeronautics, and Dr. 
Harry Bateman, Professor of Mathematics, V 
California Institute of Technology, for their 
valuable suggestions in the carrying out of this 
study and investigations. : 


APPENDIX 
EVALUATION BY CONTOUR INTEGRATION 
Cauchy’s theorem of residues 


W(m)dm= 2zi- (Sum of residues) 
Ss 

enables us to compute the definite integral in 
(20). m is complex=a+i8. The residues are 
found by integrating around each pole in the 
positive direction. Fy may be written as 


f(m)+sin 


m(1+m?)(Sinh? ym—m?) m(1+m?)(Sinh? ym —m?) 


where the meaning of f(m) is evident by com- 
paring with Eq. (20). 
Consider the complex integrals: 


M,= The imaginary part of 


emt. f(m) Adm 


M,= The real part of 


(21) 


Adm (22) 
s 


taken along the contours S consisting of the 
real axis and the infinite half circle above or 
below the real axis, see Fig. 7. 

If ¢=iog (r/a) <0 the half circle below the 
real axis must be chosen because in this case the 
integral along this part of the path converges to 
zero as m—>. For the same reason the upper 
infinite half circle is chosen when £>0, i.e., r>a. 


Consider for example the lower half-plane: Let 
m=a-—is where @ is real and positive, then, if 
0>£=—c where ¢ is real and positive, 


cos ac—i sin Bc 


= ei(a-iB)(—e) =e ica, 


eh 
which becomes zero in the limit. 


ip 
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The singularities of the integrand are all 
simple poles and are the roots of: 


Sinh? ym —m?=0 (23) 
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or if m is pure imaginary 78 
B= +sin (78). (24) 


The complex roots may be found with very good 
approximation (except m= 2) by: 


1 m(2n+1)—e (2n+1)—e 
+ +7 » 
2y 
where 


e=4 log [(2n+1)/y]/(2n+1) 


and #=2, 3, 4,5.... 
values of n: 


y=3n/2. For large 


2 4n+2 2n+1 
m= + — log +1 (27) 
3a 3 3 


The lower roots are: 


a: 0 0 0 0 +0.231 +0.320 


B: 0 +0.539 +0.909 +1 +1.629 +2.303 


where 1.629 corresponds to m=2 in (25) and 
2.303 corresponds to n= 2. 

It will be seen that for |m|=1 the roots are 
located on the imaginary axis, i.e., pure imag- 
inary. For |m| >1 all the roots are complex and 
located near the imaginary axis. 

The convergence of the power series is rapid 
and even for values of r near a, only a few of the 
poles need be considered outside the unit circle 
|m|=1. For very large || only the poles 
inside the circle |m|=1 contribute. 
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Characteristics and Functions of Thyratrons 


A. W. Hutt, Research Laboratory, General Electric Company 
(Received January 3, 1933) 


The Thyratron is a tube of very low resistance (arc 
discharge) which can be started or prevented from starting 
by a grid. Its qualities are: enormous power-amplification, 
approximately 10" per tube; efficiency between 95 and 99 
percent at all voltages above 250; unlimited size, as regards 
current-capacity; high-voltage limit equal to that of the 
Pliotron; starting time one to six microseconds; deioniza- 
tion time 10 to 500 microseconds. 

Use of Thyratron as switch and for power control. The 
controlling element may be a switch, clock, thermostat, or 
photo-tube; the controlled element a motor, magnet, 
contactor, or reactor. Typical applications of this kind of 
use are; turning on lamps at dusk, dispatching products to 
predetermined stations, cutting hot steel bars to exact 
length, opening doors at the approach of a person, wrapping 
packages, sorting beans and other articles, counting people 
or products, operating line or spot welding machines. By 


varying the phase of the grid-voltage with respect to that 
of the anode voltage, a smooth variation of average anode 
current may be obtained. A typical application is the 
dimming and blending of lights in theaters. 

Use of Thyratron as rectifier and as inverter from direct 
to alternating current. Immediate objectives are frequency 
changing, from 60 to 25 cycles for railway and power 
purposes, and from 60 to 200 for spinning mills; and 
commutatorless motors of variable and controllable speed. 
A more remote but important application is d.c. trans- 
mission of power, with inversion to a.c. at the point of 
utilization. 

Use of Thyratrons in scientific research. Applications 
include high-speed stroboscopes, timing devices, synchro- 
nous switches, and voltage regulators. Most interesting and 
promising of all is a device for counting alpha-particles, 
protons, and neutrons. 


HE Pliotron tube, or electron-tube, as it 

has come to be called, has reached the 
stage of familiarity that leaves us leisure to wish 
for something more. “If only we had a tube 
which combines the control characteristics of 
the Pliotron tube with low-resistance and large 
power capacity,” one frequently hears. 


THE IDEAL POWER TUBE 


But do you want such a tube? Consider how 
you could use it. You would connect it in series 
with a resistance or load device of some kind 
(Fig. 1), and a source of voltage. Assume it to 
be operating, carrying full load—let it be a 
typical railway load of 1000 amperes at 500 
volts. You now make the grid negative and 
reduce the current to half. The voltage drop 
across the load, which was formerly 500, is now 
250. Where are the other 250 volts? Obviously 
in the tube, whose resistance must now be high, 
high enough to absorb half the voltage and half 
the power. The grid has reduced the tube 
current by increasing the tube resistance. There 
appears to be no other way, for grid modulation 
of current and low tube resistance are mutually 
inconsistent. The difficulty is as fundamental 
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as that of the child who finds he cannot have 
his cake and eat it. We may have either low 
resistance or continuous control, by a grid, but 
not both. In large tubes, where power waste is 
important, the choice must be low resistance. 
You would not and could not use continuous 
control by a grid if you had it. 

The ideal power valve is therefore one which 
has low resistance and can be started or stopped 
by means of a grid. But would you really use 
the stopping feature? Returning to Fig. 1, let 
us make the grid more negative than before and 
completely stop the 1000 amperes. If this is 
done quickly, there will be no serious dissipation 
of power in the tube. But what will become of 


—— = 1000 AMPS. 


R L 
E=500 vorrs 


Fic. 1. Typical circuit for grid-controlled power tube. 
The grid cannot control the magnitude of the current 
without introducing high resistance in the tube. 
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the energy, }Z7, stored inductively in the motor 
armature? Physicists familiar with induction 
coils know the answer. The voltage will rise 
until all this energy is stored in the capacitance 
across the tube, or until the insulation fails. It 
is an axiom among engineers that currents should 
be “broken” only when they pass through zero. 
That is, they must be allowed to die naturally, 
and simply be prevented from restarting. Hence 
you would not use the stopping feature if you 
had it. 

I think we must conclude that the character- 
istics desired for our ideal power tube are: low 
resistance; ability to start at will, by the use of 
small control-power; and ability to remain non- 
conducting after stopping. 


DESCRIPTION OF THYRATRON 


We have such a tube, with precisely these 
qualities.! It has been named Thyratron,* from 
the Greek word Otpa, meaning door—the door 
tube. The Thyratron (Fig. 2) is a grid-controlled 
arc-discharge device. Its parts are a cathode, 
either thermionic or mercury pool type; a grid; 
anode; gas-tight envelope; and an inert gas or 
vapor, usually argon or mercury, at very low 
pressure. The pressure varies, in the case of 
mercury, from 0.001 mm to 0.050 mm, the lower 
pressure being used for high voltage and rapid 
deionization, the higher when efficiency and 
large current capacity are paramount. For argon 
and helium, the pressures should be about 3 and 
10 times higher, respectively. Slightly higher 
pressures, about twice the above values, may 
be used at low voltage. 

One structural feature deserves mention. It 
has been found possible .to use thermionic 
cathodes of a special type, which require very 
little energy to keep them hot.? This is accom- 
plished by heat-shielding the cathode on the 
outside, except for a few small openings, taking 
advantage of the fact that electrons in the 
presence of an equal number of positive tons can 


*G. E. Trademark. 

1A. W. Hull and I. Langmuir, Control of an Arc Dis- 
charge by a Grid, Proc. Nat. Acad. Sci. 51, 218-25 (1929). 
A. W. Hull, Hot Cathode Thyratrons, Gen. Elec. Rev. 32, 
213-23, 390-99 (1929). 

2A. W. Hull, Gas-Filled Thermionic Tubes, Trans. 
A.1L.E.E. 47, 753-63 (1928). 


flow around corners and through narrow crevices, 


with very little resistance, whereas heat radiation 
must proceed in straight lines. Fig. 3 shows such 
a cavity cathode, open at one end, but protected 
from loss of heat over the rest of its surface by 
several concentric metal cylinders. Inside the 
cavity are 16 barium-coated vanes, with a total 
area of 250 cm’, and normal electron emission 
of 75 amperes. The heat-loss is 60 watts. Fig. 4 
shows another form of this type of cathode, 
consisting of a large number of barium-coated 
nickel disks surrounded by heat-insulating cylin- 
ders, with a series of small holes for the escape 
of electrons. This cathode requires 300 watts to 
keep it hot and has an electron emission of 
600 amperes. An open or unshielded cathode of 
the same wattage would emit only 15 amperes. 
It is to be noted that the cavity type of cathode 
is useful only in tubes containing gas, where 
space-charge does not limit the flow of electrons. 

Now let us consider the properties of this 
new tube. 


POWER AMPLIFICATION 


In the first place, it is capable of enormous 
power amplification, thousands of times greater, 
even, than the Pliotron.* For example, a typical 
small Pliotron will give, at best, an increase of 
2 milliamperes in plate current for 1 volt change 
of grid potential. A Thyratron of the same size 
gives an increase of 2 amperes instead of 2 
milliamperes. The power required to change the 
grid potential is essentially the same in the two 
cases. This is shown in Fig. 5, which gives the 
grid current at the instant of starting in an 
FG-95 Thyratron, rated at 2.5 amperes, 1000 
volts. It is seen that the maximum current at 
starting is 0.005 microampere, at any anode 
voltage between 200 and 1000. Tests have shown 
that a resistance of 1000 megohms in series 
with the grid does not interfere with reliable 
starting. Hence the input power may be taken 
as ?R=(5X10~-*)?10°=2.510-* watts, and 
the output power at 1000 volts as EJ= 1000 X 2.5 
= 2.5 kw. The ratio (2.5 10°) /(2.5 10-*)=10", 
is the power amplification. 

The control characteristic of the FG-95 Thy- 
ratron is shown in Fig. 6. Its control ratio of 
grid to anode volatge is seen to be about 1/200, 
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Fic. 2. Mercury-arc type Thyratron. The starting of the current to each anode is controlled by a grid 
close te the anode. 


which is consistent with the above energy 
calculations. In the FG-98 Thyratron, which is 
a small filament-type tube containing argon, 
the grid current required for starting is even 
smaller, viz. 0.002 microampere over its operating 
range of 50 to 220 volts, as shown in Fig. 7. 
Its power amplification has approximately the 


same value, 10", as the FG-95. With larger 
tubes, it is probable that an amplification of 
10 times this value, viz. 10’, could be obtained 
if desired. For many purposes this high amplifi- 
cation is not needed, and “hard starting” is 
considered by some to be an advantage in power 
circuits; so that power Thyratrons are generally 
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Fic. 3. Cellular type cathode. 250 sq. cm of electron 
emitting surface are enclosed in a heat shielded cavity 
with a hole at one end for the escape of electrons. The 
shielding reduces the heating energy from 1500 watts to 
60 watts. 


designed, at present, with a power amplification 
of about 10°. 

A second important property of the Thyratron 
is its high efficiency. The first efficiency factor, 
cathode efficiency, has already been discussed. 
It may be expressed in terms of the power cost of 
electrons. This, for large cathodes, is approxi- 
mately 1/2 watt per ampere of emission. In 
practical circuits, however, the current flows for 
only 1/2 to 1/12 of the time, so that the average 


Fic. 4. Alternative type of cellular cathode. The electrons 
escape through a number of small holes. 
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. Fic. 5. The grid current required at instant of starting 
.o Thyratron, rated at 2.5 amperes, 200-1000 
voits. 
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cost, or heating power per used ampere, is from 
1 to 6 watts. 

The second efficiency factor is the voltage drop 
in the arc. A typical are-drop characteristic 
(Fig. 8) consists of 3 stages: First, a short 
portion, at approximately the ionization poten- 
tial, for small currents. Second, a ‘negative 
resistance’ transition portion, in which the 
voltage drop falls from the ionizing potential to 
about the first radiating potential. This occurs 
when the concentration of metastable atoms 
becomes sufficient, so that the ions needed can 
be produced by a two-stage process, namely 
excitation by 4.9 volt electrons and subsequent 
raising to the ionization level by impact from 
5.4 volt electrons. Third, a long positive re- 
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Fic. 6. Control characteristic of FG-95 Thyratron. The 
control ratio is approximately 1/200. 
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Fic. 7. The grid current at instant of starting in FG-98 
Thyratron, rated at } ampere, 60-180 volts. 


sistance portion, representing the increasing cost 
of producing positive ions when their number 
becomes a considerable fraction of the total 
number of atoms present. This part of the 
characteristic is a true positive resistance, as 
proven by the fact that two anodes connected in 
parallel will share the current equally, for large 
currents, under all circuit conditions.’ Practical 
tube characteristics show the first two portions 
in varying degree, according to geometry and 
vapor pressure; but the positive resistance por- 
tion is always present, if the tube is fully loaded. 
This is well illustrated in Fig. 9, which shows 
the characteristics of the FG-52 Phanotron** at 
different vapor pressures. The negative resistance 


5 
VOLTAGE DROP 


CURRENT 


Fic. 8. Typical volt-ampere characteristic of mercury 
vapor Phanotron. It consists of a short negative resistance 
portion and a long positive resistance portion. 
. 


3A. W. Hull and H. D. Brown, Mercury-Arc Rectifier 
Research, Trans. A.I.E.E. 50, 744-56 (1931). 

** The general name jPhanotron (Greek gaivw, appear, 
or be jvisible) [has been given to tubes which combine a 
thermionic cathode with low-pressure inert gas. 
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portion, which occurs at currents below 5 am- 
peres, is too short to show on this scale. 

From these data it is seen that the sum of the 
cathode and arc losses is about 7 watts per 
ampere for a mercury vapor tube with ther- 
mionic cathode** under the most favorable 
conditions. It is 3 watts per ampere for a similar 
tube with caesium vapor. The corresponding 
minimum value for a mercury-arc rectifier, or 
Thyratron with mercury pool cathode, is 20 
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Fic. 9. Volt-ampere characteristics of FG-52 Phanotron 
at different temperatures. The negative resistance portion 
of the characteristics is too short to show. 


watts per ampere, with the possibility of ultimate 
reduction to 12, and still lower if caesium is used 


in place of mercury. Under the most unfavorable 


conditions, namely high tube-drop and short 
conducting period in each cycle, the power 
consumption in the tube is between 25 and 
30 watts per ampere at full load, for either the 
thermionic or pool type. 

With these values of tube loss, it is possible 
to choose circuit and tube conditions so that the 
tube efficiency lies between 95 and 99 percent, 
at all practical voltages. 
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CURRENT AND VOLTAGE LIMITATIONS 


The current capacity of the Thyratron is 
limited only by physical size, since the positive- 
resistance volt-ampere characteristic, discussed 
above, insures proper division of current between 
the parts of a given tube, however large, or 
between any number of tubes in parallel. At 
light load, the current is concentrated in a small 
part of the tube. But as soon as that part 
attains its normal current, other parts begin 
to contribute, until all share equally at full 
load. The maximum present rating of 12-anode 
mercury-arc Thyratrons is 5000 amperes con- 
tinuous and 15,000 amperes for short periods. 
Hot-cathode Thyratrons are available with 
continuous current rating of 100 amperes; and 
there appears to be no difficulty in increasing 
the capacity to 10 or even 100 times this value. 

The voltage at which Thyratrons can be 
operated is comparable with that of high vacuum 
tubes, and the limiting value is probably the 
same for both, namely the voltage at which 
autoelectronic emission takes place, independ- 
ently of gas, due to strong electric fields. The 
highest rated voltage of present commercial 
Thyratrons is 15,000 volts, which is about the 
same as for commercial Pliotrons. But experi- 
mental Thyratrons, made by Dr. K. H. Kingdon, 
have been operated in the laboratory for long 
periods at 150,000 volts, and for short periods 
at 200,000 volts. 


TimME LAG or IONIZATION AND DEIONIZATION 


The well-known time lag in the starting of a 
spark discharge in low-pressure gases is due to 
two factors: viz., the ionization of the gas path, 
and the building up of cathode emission. In the 
Thyratron the cathode emission is ready at 
hand, and only the space-ionization needs to be 
built up; hence we may expect the time to be 
much shorter. Fig. 10 shows a cathode-ray 
oscillogram of the starting time of an FG-57 
Thyratron. The upper horizontal curve shows 
the form of voltage wave used, the lower curve 
the result of applying this voltage to the anode. 
The rise in the curve shows the time of applica- 
tion of anode voltage, and its fall, the attain- 
ment of full conduction. The sine curve is a 
timing wave. It will be seen that the (otal 


starting time, i.e., the interval between the time 
when the grid becomes positive and the time at 
which the anode current attains full value, is 
approximately 2 microseconds. This time is of 
importance for many scientific applications, such 
as the timing of stroboscopes and cathode-ray 
oscillographs, the measurement of instantaneous 
currents, and the counting of electric impulses. 


is 
Micre seconds 


Fic. 10. Time required for starting of FG-57 Thyratron. 
The two horizontal curves represent voltage between anode 
and cathode, the upper being the voltage wave without 
the anode connected, the lower with the anode connected. 
The application of anode voltage is indicated by the rise 
in the curve and attainment of full current by its fall. 
The sine curve is a timing wave. 


Deionization, or diffusion of ions to the walls 
after the current stops, requires a much longer 
time. It varies from 10 to 500 microseconds 
according to geometry and vapor pressure. This 
is the time that must elapse before the grid can 
again gain control. Fig. 11 shows the positive ion 
current to the grid of an FG-33 Thyratron, as 
measured by a cathode-ray oscillograph, im- 
mediately after stopping the anode current, for 
four values of anode current. The decay curve 
consists of two exponential portions, the first, 
of about 2 microseconds duration, representing 
the building of a positive-ion sheath around the 
grid, the second the diffusion of ions into this 
sheath. 


METHODs OF OPERATION 


On and off operation 


The simplest use of a Thyratron is as a 
switch, to turn current on and off. It has four 
properties that recommend it in preference to 
an ordinary switch or contactor, viz.: sensitive- 
ness, or smallness of power needed to control 
its operation; speed; absence of wear; quietness. 
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Fic. 11. Grid current after ‘‘stopping’’ in FG-33 Thyratron for 4 values of tube current. 


Fig. 12 shows the fundamental circuit for this 
type of service, when a.c. power is available. 
The load may be in either the a.c. or the d.c. 
branch, as indicated, and the controlling element 
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Fic. 12. Typical circuit for using Thyratron as a switch 
with a.c. power. Note that the load may be placed in 
either the d.c. (rectified) or a.c. circuit. 


may be a manually operated switch, a clock, a 
thermostat, or a photoelectric tube. The control 
voltage, E,, is shown as d.c., but a.c. is equally 
good and is more often used. The bias voltage 
E, is generally a.c. rectified by copper oxide. 


When the direct-current output is used, it 
may be employed to control motors, magnets, 
or contactors, or to saturate the d.c. windings of 
reactors in a.c. power circuits, and so control 
the current in these circuits. This type of 
operation has found many applications in in- 
dustrial control,* such as turning on lamps at 
dusk, dispatching products to predetermined 
stations, cutting hot steel bars to exact length, 
opening doors at the approach of a person, 
wrapping packages, sorting beans and other 
articles, counting people or products. 

The load may be placed in the a.c. circuit 
instead of the d.c. circuit, after closing switch A 
and opening switch B in Fig. 12. In this case, 
when the Thyratrons are not conducting, the 
transformer T constitutes a high impedance in 
the a.c. circuit; when the Thyratrons are con- 
ducting, they short-circuit the secondary of the 
transformer and make its impedance very low. 
This type of operation has found application in 


resistance welding, where its high speed and 


4W. R. King, Electron Tubes in Industry, Trans. A.1.E.E.. 
June, 1931; Chem. and Met. Eng.. June, 1932. 
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exact timing produce results hitherto unattain- 
able. 

Fig. 13 shows a method of operating a Thy- 
ratron switch with d.c. power supply. The input 
circuit and method of control is the same as 
with a.c., but in this case one Thyratron is 
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Fic. 13. Typical circuit for operating a Thyratron as a 
switch, with d.c. power. 


started to energize the circuit, the other to de- 


energize it. The second Thyratron carries only | 


a transient current, which interrupts, momen- 
tarily, the current in the first, and thus allows 
its grid to gain control. 


Phase control 


There is a still more elegant method of 
controlling Thyratrons, for which we are in- 
debted to a French engineer, P. Toulon.® It 
consists in varying the phase of the alternating 
voltage applied to the grid. This varies the time, 
in each cycle of anode voltage, at which the grid 
voltage becomes positive, thus “‘starting’’ the 
Thyratron. In the circuit of Fig. 12, if the grids 
are made positive at the beginning of each cycle 
the current will be practically continuous; but if 
they are made positive late in each cycle the 
current will flow only part of the time, viz., the 
remaining part of each cycle, and the average 
current will be less than in the first case. In this 
way one may obtain any value of average 
current, from maximum value to zero, by smooth 
continuous steps. It is to be noted that this 


5 P. Toulon, U. S. Patent 1,289,823; L. Dunoyer and P. 
Toulon, Comptes Rendus. p. 179, 1924; J. de Physique 
257-68, 289-303 (1924). 


method of controlling the magnitude of average 
current, by varying the fraction of the time 
during which current flows, is consistent with 
high efficiency. The above description applies to 
a resistance load. With an inductive load the 
current pulses will be lengthened and may 
overlap to give a continuous current, but the 
average value of this current will still depend 
on the time of starting in each cycle. 

An interesting application of this method of 
Thyratron operation is the dimming and blending 
of lights in theaters and public buildings. 


Rectification 


An important application of Thyratrons is 
controlled rectification. All practical generators of 
electric power generate alternating current, 
which must be rectified in order to give direct 
current. Thyratrons are superior to commutators, 
for this purpose, not only in their freedom from 
arcing and wear, but in their ability to control 
the magnitude of the current, thus functioning 
as voltage regulator, load limiter, and circuit- 
breaker, as well as switch and commutator; 
with the additional feature that commutation 
need not be performed at the generator, but can 
be applied at any point of the alternating- 
current network. It is believed that Thyratron 
rectifiers, thus used, have an important appli- 
cation to the problem of feeding power from the 
a.c. to the d.c. systems in large cities. 


Inversion 


The most fascinating application of Thy- 
ratrons is changing direct current into alter- 
nating. This operation is simply the inverse of 
rectification, and any good Thyratron rectifier 
can be operated in either direction by properly 
controlling the grids. Thus a transformer for 
alternating current can be converted into a 
direct-current transformer by the addition of 
Thyratrons. It is possible that electric power 
will be transmitted in the future as direct 
current, especially in cables, on long distance 
lines, or near telephone wires; or in large systems, 
where synchronization of alternating machines 
is difficult. But such large-scale applications 
must await thorough proving of apparatus and 
methods, and the working out of many details. 
One of these details is the effect of the load on 
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the inverter. Fig. 14 shows an inverter developed 
by C. H. Willis,’ which is capable of supplying 
alternating current to a load of any power 
factor, either lagging or leading. . 


Fic. 14. Inverter circuit with harmonic commutation. 
Commutating voltage is introduced by the transformer 
shown at the bottom of the diagram. (C. H. Willis). 


Frequency changing 

Frequency changing from a higher to a lower 
frequency is a special case of rectification, in 
which the Thyratron grids allow current to flow 
for a number of cycles first in one direction, and 
then in the other. Similarly, changing from lower 
to higher frequency is simply inversion, in which 
the power is drawn alternately from the three 
phases of the lower frequency. Examples of the 
first type are the change from 60 to 25 cycles 
for railway and other power purposes; of the 
second type the change from 60 to 180 for 
spinning mills. 


Commutation of motors 


The commutation of motors, like generators, 
can be performed much more perfectly by 
Thyratrons than by brushes. A special advantage 


°C. H. Willis, Harmonic Commutation Using Thyratron 


Tubes, Gen. Elec. Rev. 35,_632-8 (1932) A.I.E.E. paper © 


read January 24, 1933. 
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of Thyratron commutation is speed control. 
Fig. 15 shows a synchronous type motor with 
Thyratron commutators,’ capable of operating 
at any speed from zero up to synchronism. 


SCIENTIFIC APPLICATIONS 


The smallest application of Thyratrons, in 
physical size, may become the greatest if it 
should lead to an important discovery. This is 
their use for scientific investigation. Applications 
of this kind are just beginning, but they include 
a stroboscope® faster than any hitherto devised; 
a linear timer® for the cathode-ray oscillograph; 
a synchronous switch” for interrupting filament 
heating current while measurements are being 
made; thermostats; voltage regulators; and 
circuits for measuring short time intervals," 
currents of short duration, and end-points in 
chemical reactions. 

The most interesting scientific application is 
the counting of protons and other products of 
nuclear disintegration or cosmic radiation. Thus 
de Bruyne and Webster,” and Sparks and Picker- 


‘ for 


Fic. 15. Variable speed synchronous motor with Thyratrons 
as commutators. (C. H. Willis). 


7C. H. Willis, The Thyratron-Commutor Motor, Gen. 
Elec. Rev. 36, February, 1933. 

8H. E. Edgerton, Stroboscope Using Thyratron Tubes 
(to be published soon). 

* A. L. Samuel, Rev. Sci. Inst. 2, 532 (1931). 

1 W. B. Nottingham, Phys. Rev. 37, 1690A (1931). 

 K, J. Sixtus and L. Tonks, Phys. Rev. 37, 930 (1931). 

12 N. A. de Bruyne and H. C. Webster, Proc. Camb. Phil. 
Soc. 27, 113 (1931). 
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Fic. 16. “Scale of Two’’ Thyratron Counter, by C. E. Wynn-Williams. 


ing have applied the Thyratron to the recording 
of coincidences of two or more Geiger counters. 
Rutherford and Wynn-Williams'* have used the 
combination of vacuum tube amplifier and 
Thyratron to count alpha-particles, making 
possible the separation of particles of different 
range, or different ionizing power at any part of 
their range, to a degree hitherto unattainable. 
Most important of all, Wynn-Williams has 
devised a method (Fig. 16) of high-speed 
counting by means of Thyratrons, which appears 


%B. Sparks and W. Pickering, Meeting of American 
Physical Society, December 17, 1932. 

4 Rutherford, Wynn-Williams and Lewis, Proc. Roy. 
Soc. A133, 351 (1931). 

%C. E. Wynn-Williams, Proc. Roy. Soc. A132, 295-310 
(193i); A136, 312-24 (1932). 


capable of counting, infallibly, any number of 
particles, at a speed limited only by the de- 
ionization time of the Thyratrons. This time 
may be made less than 10 microseconds, so that 
it should be possible to register separately 
particles which arrive 10~° sec. apart, and 
perhaps as close as 10~® sec. apart. 

These methods, with better Pliotron amplifiers, 
of which the PJ-11" is a first example, and 
faster Thyratrons, of which the FP-95 is a 
forerunner (See characteristics, Figs. 5 and 6) 
will help to speed the attack on the nucleus, 
which many laboratories, with million-volt mis- 
siles, are preparing to storm. 


1% G. F. Metcalf and T. M. Dickinson, Physics 3, 11-17 
(1932). 
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Mutual Impedance of Long Grounded Wires When the Conductivity of the Earth 
Varies Exponentially with Depth 


VOLUME 4 


Marion C. Gray, American Telephone and Telegraph Company 


This paper presents a formula for the mutual impedance 
of long grounded wires above the surface of the earth, 
on the assumption that the conductivity of the earth varies 
exponentially with depth according to the formula 
y=yoe*. For b=0 the formula reduces to the known 
result for a uniformly conducting earth, while if } is 
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allowed to become negatively infinite it reduces to the 
result for an earth consisting of a conducting layer at the 
surface only. For small values of } the first terms in the 
expansion of the impedance formula in powers of 6b are 
obtained, and curves are included of the real and imaginary 
parts of the coefficient of 5. 


RMULAE for the mutual impedance of 
overhead wires parallel to the surface of the 
earth have been published with various assump- 
tions as to the conductivity of the earth. Thus 
Pollaczek' and Carson? assumed a_ uniform 
conductivity, while Mayr*® assumed a conducting 
layer at the surface only. In both these cases 
the conductivity is defined by a single parameter, 
but Haberland‘ combined the two assumptions 
to consider the two-parameter case of an earth 
of uniform conductivity with a thin surface layer 
of greater conductivity. More recently Evans® 
has generalized the Haberland case by assuming 
that the upper conducting layer is of finite 
depth, thus introducing three conductivity pa- 
rameters. In the present paper we consider 
another two-parameter case and derive a mutual 
impedance formula on the assumption that the 
conductivity of the earth varies exponentially 
with depth.® 
We consider an infinite wire of negligible 
diameter at height 4 above the surface of the 
earth (assumed plane), and carrying a current 
Te*t. As usual the exponential factor is omitted 
throughout. The surface of the earth is taken 
as the x—y plane, the wire being parallel to the 


1F, Pollaczek, Elekt. Nachr. Tech. 3, 339-359 (1926). 

2 J. R. Carson, Bell Sys. Tech. J. 5, 539-554 (1926). 

*0. Mayr, Elektrotech. Zeits. 46, 1352-1355 (1925). 

*G. Haberland, Zeits. angew. Math. 6, 366-379 (1926). 

5H. P. Evans, Phys. Rev. 36, 1579-1588 (1930). 

* A similar assumption has been made recently by G. J. 
Elias, Elekt. Nachr. Tech. 8, 4-22 (1931), in connection 
with a different problem. 
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x axis, and the z axis extends vertically upwards 
into the air. The frequency, w/(27), is assumed 
to be sufficiently low to allow all displacement 
currents to be neglected. The inductivities of 
both earth and air are taken equal to that of 
free space v, and their capacitivities zero; the 
conductivity of the air is also assumed to be 
zero, while the conductivity of the earth is 
defined by the formula 


v= 20, 


(1) 


where yo is the conductivity at the surface 
(z=0), and we shall at first assume that the real 
constant 0 is positive. All distances are measured 
in meters, impedances in ohms and conductivities 
in mhos per meter; v has the value 1.256 10~° 
henries per meter. 

Since the wire is very long we may neglect 
the components E, and E, of the electric field 
as compared with £,, and the differential 
equations satisfied by E, and the magnetic 
components //, and J/, are therefore, 


iwvll,= —0E,/02, (2) 
iwvll,= dE,/dy, (3) 
E,= (4) 


where the electric field is expressed in volts per 
meter, the magnetic field in ampere-turns per 
meter. From these equations the differential 
equation for E, is easily found to be 


(5) 
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where I'o= (iwvyo)'” is the propagation constant 
at the earth’s surface. A solution of this equation 
which vanishes at an infinite depth is 


E.= cos yu du, (6) 


where we have written 8 for 1/b. Here K.(¢) is 
the Bessel function of the second kind with 
imaginary argument as defined by G. N. Wat- 
son,’ and f(z) is an arbitrary function of its 
argument to be determined by the boundary 
conditions. These conditions are that E, and H, 
must be continuous at the surface, z=0, and 
from Eq. (2) it follows that the continuity of /7, 
may be replaced by the continuity of dE,/dz. 

In the air the electric field E,9 at the point 
(y, 2) may be formulated as 


+f yu du, (7) 


where the first term gives the field due to the 


primary current in the wire,’ while the second 
term represents the secondary field, ¢(u) being 
at present arbitrary. From Eqs. (6) and (7) 
using the two boundary conditions, we have the 
following equations to determine the two func- 
tions f(u) and ¢(u), 


0)f(u) = o(u), (8) 
ToK 0) f(u) = (tw (9) 


where in obtaining Eq. (9) we have used the 
Fourier integral 


h/(l?+y*)= e~™ cos yu du. 
0 


With the use of the recurrence formulae for the 
Bessel functions, the solution of these equations 
is 


0) 
Kogy(2BTo) 


= — (11) 
0) 


f(u)= 


(10) 


The electric field in the air is therefore determined as 


Ew=— 


e~ Cos yu du. (12) 
K 0) 


The mutual impedance per unit length between two infinite parallel wires at heights h and JT 
above the earth’s surface may be written down immediately: 


Qe 


Kog,(28T 0) 
f = e~ cos yu dp. (13a) 
0 


K 0) 


The solution (13a) has been obtained on the assumption that the conductivity of the earth is 
increasing exponentially with depth (0<b). If, however, we wish to include negative values of b, 
the formula is easily modified, the only change in the derivation being that the solution (6) of Eq. 
(5) is no longer valid, and must be replaced by the solution 


E.= f(u) I-26 u( — cos yp du, (6’) 


where J,(¢) is the Bessel function of the first kind with imaginary argument as defined by 
Watson.’ Then in place of (13a) we find the formula 


Z2=— 


og 
Po 


7G. N. Watson, Theory of Bessel Functions, Cambridge, 


1922, page 78. 


—26T») 
= = Cos yp dy. (13b) 
0 


— 2610) 


8 Reference 5, page 1583. 
® Reference 7, page 77. 
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If we put b=0 both (13a) and (13b) reduce to the well-known formula” 
log 


for a uniformly conducting earth. Further, by allowing 8 to approach zero through negative values 


2)12— cos yu du (14) 


while at the same time yo becomes infinite in such a way that y;= — Byo remains finite, we obtain 
log += 008 yu dy, (15) 


which is Mayr’s result" for a conducting layer at the surface only. 

Returning to Eq. (13a) let us simplify it by making the substitutions $’=I9(2z)—/,, 
y’ = T'(2i)-'*y, with similar definitions of h’ and H’, thus multiplying all distances by the attenu- 
ation constant. We then have 


e cos y'u du. (13) 


For the important case of small values of b’(=1/8’) it is possible to expand Z):2 in a series of 
powers of b’, the first terms of the expansion being 


Z12= (wv/2n) log Ly? +(h! +11’)? + — Poh’ +11, y’) 
+0'Pi(h' +H’, y’)+b?Po(h' +I’, y')+-++}, (16) 


where 
Cut cos yu dy, (17) 
0 
i(h, oJ, cos yu du (18) 
Pa(h, y) 4y?+ 3% 4p (19) 


This expansion is valid for either positive or negative values of b’. In Eq. (16) the first term is the 
mutual impedance for an earth of infinite conductivity, the second term is the Carson correction 
for a uniform finite conductivity, while the third term gives an approximate correction for an earth 
whose conductivity increases or decreases slowly with depth. 

An expansion of the integral (17) in a form suitable for computation has already been given 
by Carson” in terms of the polar coordinates (r, 0), where h=r cos 0, y=rsin 6. A similar, though 
more complicated, expansion for the integral (18) is as follows: 


Ph, y) =P,(r, U\(r, 0)+iVi(r, 6) 
where 
2U,(r, 6)=(x/8)[cos A cosh B+cos B cosh A+sin A sinh B+sin B sinh A ]—S,— 
+4(S4—S3) +3(Ss+ Ss) log 2+(S;—Se) log r— Tot (20) 
1 Reference 1, formula (23b); reference 2, formula (24); reference 4, formula (18). 


" Reference 3, formula (13). To obtain Mayr’s formula we have to put 4=0 in Eq. (15). 
® Reference 2, formulae (32) and (33). 
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2V,(r, @)=(2/8&)[cos A cosh B+cos B cosh A —sin A sinh B—sin B sinh A ]—S,+S2 


+3(S3+S;) log 2+(S3—S;) log —S3')0 +} 


(21) 


In these expansions we have written A=r(cos 6+sin 0), B=r(cos 6—sin 6), and the convergent 


series are defined as follows, 


cos 40 cos 84 
S,;=1- + 
3?-5?-7?-9 
2r? cos 20 cos cos 100 
11 
cos 50 24r? cos 90 
S3=r cos meee, 
5! 9! 
cos 50 cos 96 
Y3=rcos 6 ¥(2)— 6) - -, 
5! 9! 
r cos 56 cos 96 
S,=r cos 0— + 
272!3! 24415! 
r cos 56 E cos 96 1 
2) 22213! 6} 24415! 10 


cos 30 cos 70 2°r'!! cos 110 


3! 7! 11! 


cos 30 cos 70 


3! 


r'! cos 116 


2°5!6! 


r> cos 30 cos 70 


+ 
2-1!2! 23314! 


+ 


r® cos 30 1) 1’ cos 70 
2) 
2-1!2! 23314! 
and S,/(r, 0)=S,(r, n=3, 4, 5, 6. 


The function ¥(m) appearing in the ~ series 
is the logarithmic derivate of the gamma- 
function." 

From these expansions curves have been drawn 
showing the real and imaginary parts of P(r, @) 
as a function of r for values of cos @=0, 4, 3, 3, 1. 
Figs. 1 to 4. Equal intervals of cos @ were used 


3 See e.g., E. T. Whittaker and G. N. Watson, Modern 
Analysis, Third Edition, Cambridge, 1920, page 241. 
When is an integer ¥(n +1) = —y+1+34+ 
where y is Euler’s constant and has the value 0.57722. 


2°r"! cos 1186 


1 r'! cos 110 1 

|wo+—]- 
8 25516! 12 


rather than equal intervals of @ itself, since it 
was found that this procedure gave a much more 
evenly spaced set of curves. For convenience of 
reference, curves for the real and imaginary parts 
of P(r, @) are also given, these being obtained 
from Carson’s expansions with slight alterations 
due to changes of notation. It may be worth 
while to note that for wires on the surface of 

4 Tn drawing the curves for cos @=1 and cos @ = }, use was 


made of the tables of values published by J. E. Clem, 
Trans. Amer. Inst. Elect. Eng. 50, 909-910 (1931). 
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wusor 


Fic. 1. U;(r, 6) = Real part of P,(r, @). 


“TT | | = | | 
| | | 
| | | 
AN 
Fic. 2. Vi(r, 6) =Imaginary part of P,(r, 6). 
the earth (@= 7/2), the expansions (20) and (21) nmi _— 


take the much simpler form Fic. 4. Vo(r, 6) =Imaginary part of Po(r, 4). 


y+sin y) 
+2-!?rber’(2'/2y) — Si(y, 0) +S2(y, 0), | functions of the first kind and of complex 
, Finally for large values of r an approximation 
— 
— Si(y, 0) — 0), to the value of @) is given by the asymptotic 
where ber’(z), bei’(z) are derivatives of the Bessel | expansion 


cos @ cos2@ 3.cos4@ 3cos56 
Pi(r, = + + 
4r 4r? 


cos 20 cos3@ 3cos 
4r8 8r4 


In conclusion I wish to record my gratitude to Mr. R. M. Foster for many helpful suggestions 
in the course of this work. 


% Tables of values of these functions have been published by H. B. Dwight, Trans. Amer. Inst. Elect. Eng. 48, 812 
(1929). 
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Fic. 3. Uo(r, 0) = Real part of Po(r, 0). 
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Capacitance and Potential Gradients of Eccentric Cylindrical Condensers 


CuesteR L. Dawes, The Harvard Engineering School 
(Received November 7, 1932) 


OMETIMES it becomes desirable to calcu- 
late the capacitance and the potential 
gradient of eccentric cylindrical condensers, for 
example, with cables and with special condensers. 
Although methods have been evolved which give 
the capacitance and potential gradient in such 
condensers *:42 very little literature giving a 
complete but compact treatment of the subject is 
available. Moreover, the different formulas pre- 
sented here are derived on a somewhat different 
basis than is ordinarily adopted. Several forms of 
the same formula together with their analysis are 
presented. The formulas are derived from the well- 
known method of inverse points and images.’ For 
simplicity, the formulas are derived in statfarads 
per cm length, with vacuum as the dielectric 
(kx=1). It is a simple matter to apply the 
necessary coefficients for the desired lengths and 
for other dielectrics. 

Fig. 1 shows the cross section of the two 
eccentric cylinders, both of which are very thin; 
the radius of the outer cylinder is 6 cm and its 
center is at C’; the radius of the inner cylinder 
is a cm and its center is at C. The distance 
between C’ and C is c cm. Let P be any point 
on the surface of the inner cylinder; draw CP. 


Fic. 1. Eccentric cylinders and inverse points. 


* Reference 3, pp. 168-170. 

‘Mayo D. Hersey, The Resistance, Inductance and 
Capacity of Eccentric Cylinders. Elec. World 56, 435 (1910). 

2 Alexander Russell, Problems in Connection with the 
Electrified Cylindrical Conductors. J. Inst. Elec. Eng. 64, 238 
(1925-26). 

3 Alexander Russell, A Treatise on the Theory of Alter- 
nating Currents. Cambridge University Press (1914). 
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Also draw the straight line BCC’ and locate 
points A and B such that 


CA XBC=a? (1) 


(2) 


The charge on the inner cylinder is +Q e.s.u. 
per cm length and the charge on the outer 
cylinder is —Q e.s.u. per cm length. Let the 
potential of the inner cylinder be V; statvolts 
above earth potential and the potential of the 
outer cylinder V2 statvolts above earth po- 
tential. 

A and B are the inverse points of the system. 
So far as the outer cylinder is concerned the 
effect of the charge on the inner cylinder is as 
if the charge +Q were concentrated uniformly 
along the axis corresponding to the point A; so 
far as the inner cylinder is concerned the effect 
of the charge on the outer cylinder is as if the 
charge —Q were concentrated uniformly along 
the axis corresponding to the point B. This may 
be demonstrated as follows: 

Let O, the midpoint between A and B, be at 
zero or earth potential. Draw OP; AP=17; and 
BP=r.. The potential at P, or the work done 
in carrying a unit charge from O to P against 
the forces due to +Q at A and —Qat B, is 


OA 2 T2 20 
W=V,= —dr+i —dr 


" r OB r 
= 20 log. (3) 


where log, is the logarithm to the Napierian base. 
In triangles CPA and CPB the angle @ is 
common from (1), CA/CP=CP/BC; they are, 
therefore, similar. 
It also follows that 


BP/AP=r2/r,= BC/CP. 


and 


(4) 


Maintain r2/r; constant so that V, remains 
constant. 
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It is well known that, if a point so moves that 
the ratio of its distances from two fixed points is 
constant, its locus is a circle. Hence the locus of 
P is a circle, since from (1) and (4) CP or its 
radius a is constant. From Green’s theorem, the 
presence of the cylinder with a as its radius does 
not change the potential relations due to the 
charge along the axis corresponding to the point 
A, at points external to the inner cylinder. 

Let P’ be a point on the inner surface of the 
outer cylinder; draw C’P’, AP’ and BP’ where 
AP’=r,' and BP’=r,’. Triangles BC’P’ and 
AC’P’ have the angle 6’ in common and from 
(2), C’A/C’P’=C’'P’/BC’; these triangles are, 
therefore, similar and 


BP’ /r;'= BC'/C'P’. (5) 


If ro’/r;’ is maintained constant, the locus of 
P’ is also a circle with C’ as a center and } as 
a radius. In a manner similar to that used to 
determine the potential of point P (3), it may 
be shown that the potential of P’, 


V2=2Q log. (12’/r1’). (6) 


By Green’s theorem the presence of the outer 
cylinder does not disturb the potential relations 
in the field due to the charges at A and B. 
Therefore, so far as the inner cylinder is con- 
cerned, the charge on the outer cylinder acts as 
if it were concentrated along the axis corre- 
sponding to the point B. 

Let L and M be the intersections of the line 
AB with the inner and outer cylinders respec- 
tively; let the distance CA=p and the distance 
C’A=p’. 

When point P coincides with Z and point P’ 
with M, from (4) and (5), 


BC/a; /r;'= BC'/b. 
From (3) and (6), 
Vi=2Q log. (BC/a); V2=2Q log. (BC’/b) 
Vi — 20 log. (b/a) X (BC/BC’). 


C= 


Vi-Vse 


m—(m?—4a%2)"2 


From (1) and (2), 
BC=a*?/CA=a?/p; 


Hence, 


Vi— V2=2Q log. ((b/a)(a*p’/b*p)) 
and the capacitance 
C= Q/(Vi-— V2)=1/[2 log. ((b/a)(a*p’ /b*p)) (7) 


Before the capacitance can be found, p and p’ 
must be determined. 


BC=a"/p; BC’=B*/p’ 
BC’—BC=c B/p’—a*/p=c (I) 
c= p'—p. (ID) 
Solving (I) and (II) simultaneously, 
p=m/2c+[(m/2c)?—a?}'”, (8) 
p’=m' (9) 


where m= 0?—a?—c? and m’=8?—a?+c. 

When the minus signs in (8) and (9) are used, 
the values of p and p’ are equal to CA and C’A. 
If the positive signs are used, the values of p 
and p’ are equal to CB and C’B, that is the 
distances from C and C’ to the other inverse 
point B. The correct result is obtained with 
either sign; with the positive signs, the value of 
capacitance becomes negative, corresponding to 
a direction of work which is negative to that 
assumed in deriving (3) and (6). 

The capacitance may be determined by sub- 
stituting the values of » and p’ of (8) and (9) 
in (7). 

The capacitance may be expressed in four 
quite different forms. The form which may be 
used is determined by whether it is convenient 
to use logarithmic functions or hyperbolic func- 
tions or whether it is desired to make comparison 
with the capacitance of a concentric cylindrical 
condenser in which the cylinders have the same 
radii as for the eccentric case. 

(a) By substituting (8) and (9) directly in (7), 


This formula involves logarithmic functions only. Positive signs before the radical may be used as 


well as negative signs. 


(b) The substitution of (8) and (9) in (7) may be made as follows: 
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b a 
C= {2 log. 2 log. 
a 


(m!? — 2 4) -1 
m' —(m c?) | a1) 


m— (m? 4a*c?)"/2 


This formula is advantageous in that the 
capacitance of the system when the cylinders 
are concentric is [2 log, (b/a)]}~“ statfarads per 
cm length. Hence the second term of the ex- 
pression is a correction term for the eccentricity 


of the cylinders. Thus the effect of eccentricity 
on the capacitance of a cylindrical condenser is 
directly determined. 

(c) The substitution in (7) may also be made 
as follows: 


b a® b(m'/2be+((m' 4 


a P 


a(m/2ac+((m/2ac)?—1)"? 


= {2[log. (m’/2be+((m' /2bc)? —1)"*) —log. (m/2ac+ ((m/2ac)? 
= {2[cosh—(m’/2bc) —cosh~! (m/2ac) (12) 


Hence the capacitance of eccentric cylindrical 
condensers varies inversely as the difference of 
two hyperbolic angles. 

Because of the use of the positive signs in (8) 
and (9), the value of the capacitance determined 
by (12) is negative. 

A. E. Kennelly* has derived the equation for 
the capacitance in a form which is quite similar 
to (12). 

As ¢ approaches zero, each of the two hyper- 
bolic angles in (12) approaches infinity, but their 


difference always remains finite. With decreasing 
values of c, this difference increases at a di- 
minishing rate, approaching as a limit a constant 
value corresponding to the capacitance of the 
cylinders when they become concentric. 
(d) Eq. (12) may be further simplified. From 
(8) and (9), 


ra (1/2c)[m — (m? —4a?c?)"/? 
By expanding, m” —4b*c?= m?—4a°c?. Hence, 


Expanding m and m’, 


2a? 
Substituting in (7), 
C= {2 log. 
where B= (b?+a?—c?)/2ab. Hence, 
C=[2 B}". (14) 
In the foregoing analysis, the capacitance has 


‘A. E. Kennelly, The Linear Resistance between Parallel 
Conducting Cylinders in a Medium of Uniform Conductivity. 
Proc. Am. Phil. Soc. XLVIII (1909). 

5 Alexander Russell gives Eq. (13). Reference 1, 1904 
edition, page 104. 


p’ (b?+-a?—c?) + [(b? +a? —c?)? b | 


a 2ab 2ab 


been given in statfarads per cm length with 
vacuum as the dielectric. With a dielectric whose 
permittivity is x, the microfarads per 1000 ft. 
are obtained by multiplying by 3.388-10-°-«; 
likewise microfarads per mile are obtained by 
multiplying by 17.89-10-*-«. 

As an example of the application of (14), let 
b=2 cm; a=0.5 cm; c= 1.0 cm; and «x= 1.0. 

The calculated capacitance is 0.5686 statfarad 
per cm length and the capacitance, when the 
cylinders are concentric, is 0.3606 statfarad per 


cm length. The ratio of the capacitance of this 
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concentric arrangement is 1.30. 


THE POTENTIAL GRADIENT 


Pender and Osborne® have derived the fol- 
lowing Eq. (15) giving the potential gradient at 
the surface of either of the cylinders of a con- 
denser consisting of two equal parallel cylinders. 
The cross section of such a condenser is shown 
in Fig. 2. The radius of each cylinder is a cm. 
The distance CC’ between the axes of the 
cylinders is D cm, and the distance between the 
inverse points A and B is d cm. 


Fic. 2. Parallel equal cylinders and inverse points. 


The potential gradient at the surface of either 
cylinder at P, the intersection of any radius r, 
is given by 


dV /dr=2Qd/[xa(D—2a cos 6) (15) 


where Q is the charge per cm length, «x is the 
dielectric constant, and @ is the angle which the 
radius r makes with line CC’. Eq. (15) is the 
polar equation of an ellipse with the origin at 
one of the foci. 

So far as any effect on cylinder (I) (Fig. 2) is 
concerned, the charge on (II) may be assumed 
as being concentrated uniformly along the axis 
corresponding to the inverse point B. This 
system becomes analogous to that of Fig. 1, 
since so far as any effect on the inner cylinder 
is concerned, the charge on the outer cylinder 
may be assumed to be concentrated along the 
axis corresponding to the inverse point B. 
Hence, (15) is applicable to the system shown 
in Fig. 1. It only remains to convert the param- 
eters Q, d and D into those corresponding to 
the system shown in Fig. 1. 


® Harold Pender and Harold S. Osborne, The Electrostatic 
Capacity between Equal Parallel Wires. Elec. World 56, 667 
(1910). 
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particular eccentric arrangement to that of the 


Referring to Figs. 1 and 2, from (1) 
BCX CA=a’, 
(d+p)p=a’, 
where p in Fig. 2 equals CA 
d= (a’—p*)/p (16) 


where d= AB in Fig. 
In Fig. 2, D=d+2p. Hence, 


D= (a?+p")/p. (17) 
Substituting (16) and (17) in (15), 
V 
(18) 


Or « a(a?+p?—2ap cos 
Q=CE, where E=V,-Vs2. 


C is most conveniently given by (14), although 
any of the other Eqs. (10), (11) and (12) may 
equally well be used. 


Hence 
KE 
Q= = (19) 
2 cosh! 8 2cosh—[(b?+a?—c?)/2ab] 
and 
fa] V 
—=E - (20) 
or a(cosh- 8) (a?+ p?—2ap cos @) 
=E (21) 


ar,?(cosh™ B) 


Hence the potential gradient varies inversely 
as the square of the distance of the point P from 
the inverse point A. 

For example, if a=0.5 cm; b=2 cm; c=1 cm; 
and E=10,000 volts, p=0.094 cm, and from 
(20), 


dV 10,000(0.25 —0.0088) 


0.5(cosh™ 1.625) (0.25 +0.0088 — 0.094 cos 4) 
volts per cm 


= 48,100/(2.755—cos @) volts per cm. (22) 


The maximum gradient occurs when 6=0 and 
is equal to 27.4 kv per cm; the minimum gradient 
occurs when @= 180° and is equal to 12.8 kv per 
cm. The corresponding value of the potential 
gradient with concentric cylinders is 14.42 kv 
per cm. 


CAPACITANCE AND POTENTIAL 


Fic. 3. Polar diagram of potential gradient. 


Fig. 3 shows a polar graph of the potential 
gradient represented by (22). 


THE POTENTIAL GRADIENT AT ANY POINT 
WITHIN THE DIELECTRIC 


Let it be desired to determine the potential 
gradient at any point P; within the dielectric of 
an eccentric cylindrical condenser, Fig. 4. As in 
Fig. 1, C is the center of the inner cylinder 
whose radius is a. C’ is the center of the outer 
cylinder whose radius is 6 and A and B are the 
inverse points of the system. 

Draw P,C=r, P,A and P,B; let @ be the 
angle between P,C=r and CB. If P, is deter- 
mined by the polar parameters r and 6, P\A 
and P,B are readily determined trigonometri- 
cally, since CA=p and AB=d are already 
known. Draw such that <C,P,A= <P,BA, 
C, lying on BA extended. C; is the center of 
the equipotential circle passing through P; and 
C\P, is its radius a;. Let M, be the intersection 
of the circle with AB; also let C,A= py. 

Triangles C;P,A and C,P;B are similar, and 


C\A [C,M,= pi/ai= C,\M,/C,B, 
(23) 
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Fic. 4. Potential gradient within dielectric of eccentric 
cylindrical condenser. 


It also follows that 
AP, /BP,. (24) 
Solving (23) and (24) simultaneously, 
(25) 
a,= (26) 


The angle @, between the radius P,C; and 
is readily determined since P;A, 
and C,A= ;, are all known. 

Substituting (19) in (18), making a=a,, p= pi 
and 6= 6; corresponding to point ~;, 


OV a’—p; 
(27) 
dr’ a,(cosh™ cos 41) 


a ini” [cosh~! 


where AP). 


For example, under the conditions from which 
(22) was obtained, with r=1.4 cm, 6=135°, 
computation gives AP,=r,;'’=1.467 cm; BP, 
= 3.785 cm; p:= 0.4537 cm; a,=1.171 cm; and 
122.6°. 

dV /dr=4.335 kv per cm as compared with 
13.89 kv per cm at the surface of the conductor 
at its intersection with the same radius 7. 

The author is indebted to Dean H. E. Clifford 
of The Harvard Engineering School for his many 
helpful suggestions in the preparation of this 


paper. 
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